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ABSTRACT 

Extensions of Einstein gravity with quadratic curvature terms in the action arise in 
most effective theories of quantised gravity, including string theory. This article explores 
the set of static, spherically symmetric and asymptotically flat solutions of this class of the¬ 
ories. An important element in the analysis is the careful treatment of a Lichnerowicz-type 
‘no-hair’ theorem. From a Frobenius analysis of the asymptotic small-radius behaviour, 
the solution space is found to split into three asymptotic families, one of which contains 
the classic Schwarzschild solution. These three families are carefully analysed to determine 
the corresponding numbers of free parameters in each. One solution family is capable of 
arising from coupling to a distributional shell of matter near the origin; this family can 
then match on to an asymptotically flat solution at spatial infinity without encountering 
a horizon. Another family, with horizons, contains the Schwarzschild solution but includes 
also non-Schwarzschild black holes. The third family of solutions obtained from the Frobe¬ 
nius analysis is nonsingular and corresponds to ‘vacuum’ solutions. In addition to the three 
families identified from near-origin behaviour, there are solutions that may be identified 
as ‘wormholes’, which can match symmetrically on to another sheet of spacetime at finite 
radius. 
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1 Introduction: Second plus fourth-order gravity 


The inclusion of quadratic curvature terms into the gravitational action is principally moti¬ 
vated by the form of one-loop quantum corrections [lj. In 4D spacetime there are effectively 
only two independent quadratic-curvature integrated invariants, owing to the existence of 
the Gauss-Bonnet topological invariant. Starting from the correspondingly general second- 
plus-fourth-order action^] 


I = d 4 xyf^ (jR - aC^ prT C^ pa + PR 2 ) 


( 1 . 1 ) 


(in which C pU p a is the Weyl tensor, i.e. the traceless part of the curvature tensor R pU pa), one 
obtains^] a renormalisable system 12 . The spectrum of this theory contains [3] a massless 
graviton, a massive spin-two ghost excitation with (m 2) 2 = and a massive non-ghost 

= where G is 


spin-zero excitation with (mo ) 2 = Pa ■ The canonical value of 7 is 


6/3' 


1 

I 67 rG 


the 4D Newton constant. 


The renormalisable quantum system (1.1) is also asymptotically free mu in the sense 
that if one writes the coefficients of the quadratic-curvature terms in Yang-Mills style as 
1 /<72 and 1 /, then both couplings 52 and go tend to zero at large energies. This raises 


the question as to whether the high-energy regime of the model ( 1 . 1 ) might avoid the 


problems associated with the spin-two ghost in the spectrum by effectively decoupling that 
excitation at high energies. Such issues have recently been discussed in the context of the 
asymptotic safety program for quantum gravity [ 6 ], but to date there does not appear to 
be a consensus on this point. A key problem in this approach is to obtain robust results 
that are not renormalization-scheme dependent. A related question, already at the classical 
level, is whether the interaction structure of the theory might even be such as to avoid the 
destabilisation of the vacuum by ghost-driven instabilities (7]. 

Gravitational theories including quadratic curvature terms arise generically in all ap¬ 
proaches to quantum gravity. In particular, the Gauss-Bonnet combination 


Igb = f V^g(Rp» P aR^ pa - ±RpvR pv + R 2 ) 


( 1 . 2 ) 


Using the 4D Gauss-Bonnet theorem (1.21, the action may also be written as 


J d' l y/^(' r R-2aR^R^ + (/3+^)R 2 ). 


Strictly speaking, for renormalisability one should also include a cosmological constant in (1.11. Note that 


the parametrisation of the higher-derivative terms in the action differs from that used in Refs EiE]. 
Specifically, o^gj-g = 2 a Ref. (2j an< ^ ^here = ^Ref. |2] ~ 3 Q Ref. [2]' 
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is a topological invariant in four spacetime dimensions, but not in higher dimensions, where 
it falls into the class of Lovelock terms [SJ. It occurs in the quantum effective action of 
heterotic string theory in 10 spacetime dimensions |9]. Various styles of dimensional com- 


pactification of (1.2) can then yield the quadratic terms of (1.1) in a variety of combina¬ 
tions]^] In dimensions D > 4, the Lovelock-Gauss-Bonnet combination (1.2) also allows for 
cosmological solutions flO) . However, in order to keep our considerations clearly focused, 
we shall restrict our attention in this paper purely to four-dimensional spacetime gravity 


derived from (1.1) without a cosmological constant. 

In this paper, we will not be concerned with difficult questions of the full physical 


acceptability of the theory (1.1) at the quantum level. Instead, we shall adopt a working 
assumption that, in whatever emerges as an acceptable quantum theory of gravity, the 


system (1.1) may be a dominant part of the effective action at least for some ultraviolet 
scale of energies. This might have, for example, cosmological implications, which could in 
turn indicate a scale for the quadratic-curvature term coefficients. It might also be the case 


that the effects of the quadratic-curvature terms in (1.1) are also characteristic of those of yet 
higher-order terms. Whatever the fate of the negative-energy massive spin-two excitation, 
we shall adopt the point of view that its effect on static classical solutions should nonetheless 


be considered. Accordingly, we shall adopt the action (1.1) as is, and shall consider the 
implications of its field equations for spherically-symmetric static solutions. We shall thus 
treat the fourth-order terms on an equal footing with the second-order terms, and not just 
as perturbations to the Einstein theory. 

Some aspects of the classical solutions to the second-plus-fourth-order gravity theory 
are well-known. In Ref. [3], an analysis was given of spherically-symmetric solutions in the 


linearised limit of the theory (1.1) when coupled to point and extended sources. As one can 
expect from a theory whose dynamical spectrum involves massive spin-two and spin-zero 
modes as well as the massless spin-two Einstein mode, the static solutions to the linearised 
theory involve both a £ potential arising from massless spin-two virtual particle exchange 
and ^— Yukawa potentials arising from m = m 2 massive spin-two and from m = mo 
massive spin-zero virtual exchanges. Moreover, by writing the spherically-symmetric and 


In the process of dimensional reduction, various massless scalar fields are generated which combine with 
the D = 10 dilaton. When the dimensionally reduced D = 4 theory is written in Einstein frame, scalar 
held prefactors appear in front of the curvature-squared terms in the effective action. The study of string- 
generated higher-derivative gravity models accordingly requires consideration of such scalars together with 
the quadratic curvature terms. In this paper, however, we shall restrict attention to purely geometric terms 
in the action. 
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static spacetime metric in Schwarzschild form 


ds 2 = —B(r) dt 2 + A{r) dr 2 + r 2 d0 2 + r 2 sin 2 9 d<p 2 , 


(1.3) 


assuming a Laurent expansion of A(r) and B(r), and carrying out a Frobenius-method 
analysis of the indicial equations for the leading asymptotic behaviour as r —> 0 in the radial 
coordinate r, it was found in Ref. [3] that the leading asymptotic behaviours]^] A(r) ~ r s and 
B{r) ~ r 4 can arise in three distinct solution families: (s,t) = (2,2), (1,-1) or (0,0). At 
the time, an initial analysis of the number of free parameters characterising these indicial 
families was made, but in a pre-computer-algebra era, the full picture of such parametric 
dependences was not easily to be found. 

In this paper, we return to a detailed study of the spherically-symmetric solutions to the 


field equations following from the action (1.1). Although the classic Schwarzschild solution 


of Einstein’s theory (which belongs to the (1,-1) family) clearly remains a solution to 


the higher-derivative theory derived from the action (1.1), we shall find that this is not a 


solution that arises from normal minimal coupling to ordinary ghost-free matter. Instead, 
we find that solutions that can arise from such ghost-free matter coupling belong to the (2, 2) 
indicial family of solutions. Subject to the additional assumption of asymptotic flatness as 
r —> oo at spatial infinity, we find that such solutions do not have a horizon, but have a 
naked singularity as r —> 0. This agrees fully, moreover, with numerical calculations of such 
solutions made in the case m 2 = mo in Ref. HU. 

If one overlooks the issue of source coupling, which in any case has been a delicate 
subject in general relativity for decades mm, then the ‘black hole’ solution family in¬ 
cluding a horizon can be investigated in its own right. Assuming in addition asymptotic 
flatness at spatial infinity, the analysis is made much simpler by a Lichnerowicz-style ’no¬ 
hair’ theorem m for the trace of the higher-derivative field equations, which implies that 
the existence of a horizon together with the assumption of asymptotic flatness leads to the 
requirement that the Ricci scalar must vanish: R = 0. Analysis of the remaining traceless 
components of the held equations is more subtle. In an earlier paper m, we reported our 
disagreement with the traceless-equation analysis of Ref. m, which would have consider¬ 
ably simplified the study of the black-hole family. In the absence of a traceless-equation 
Lichnerowicz theorem, an alternative option is to make a perturbative analysis of the black- 
hole family of solutions starting from the classic Schwarzschild solution. We obtain in this 


In context, there should be no confusion between the indicial exponents t and s here and the coordinate t 


and interval s in (1.3 1 . 
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way a result that the Schwarzschild solution is at least generally isolated , in the sense that, 
for spherically-symmetric static solutions possessing a horizon, solutions perturbatively dif¬ 
ferent from Schwarzschild necessarily must violate the condition of asymptotic flatness at 
spatial infinity. 

The general perturbative isolation of the Schwarzschild solution within the indicial 
(1,-1) solution family does not exclude the possibility of other asymptotically-flat and 
spherically-symmetric solutions with horizons that differ from Schwarzschild by a non¬ 
infinitesimal amount in the (1,-1) family parameter governing the ‘non-Schwarzschild’ 
structure of the solutions. Indeed, in Ref. m we demonstrated that this possibility is in¬ 
deed realised: there exists a range of values for the black-hole horizon radius ro, bounded 
below by a certain multiple of the I /7712 = y/2a/^ length scale, for which one obtains a sin¬ 
gle static black-hole solution in addition to the Schwarzschild solution. The corresponding 
existence of a minimum value for ro in comparison to the y/a/j scale size in the perturba¬ 
tive no-hair theorem dovetails with the numerically found existence of a branch point for 
black-hole solution phases. As one approaches this branch point, clearly the perturbative 
isolation of the Schwarzschild solution must break down. 

We begin in Section [2] with a review of the structure of the gravitational field equations 
following from the action ®> when restricted to the case of spherically-symmetric and 
static solutions, initially without considering contributions from sources. In particular, we 
discuss the reduction of the differential order of these ‘almost vacuum’ equations to get 
a better fix on the maximum number of integration-constant parameters determining a 
particular solution family. We shall find that such ‘almost vacuum’ equations reduce to a 
pair of third-order ordinary but coupled and quite nonlinear differential equations for A(r) 
and B(r). The full details of these equations are given in Appendix |A} Next, in Section 
[3j we complete the analysis of the parametric dependence of the various indicial solution 
families begun in Ref. [3j; the advent of Mathematica now makes this much more tractable. 

Given that one is principally interested in solutions that are asymptotically flat as 
r —> 00 , in which limit a linearised analysis of the solution families becomes appropriate, in 
Section [4] we next consider the spherically-symmetric static solutions to the field equations 
when linearised in A{r) and B(r). In part, this reviews the linearised solutions found already 
in Ref. [3j, but with a key addition: we now consider in some detail the matching between 
an interior vacuum and the exterior solution when matching across a shell delta-function 
source. 

Coupling to shell delta-function sources in the full nonlinear theory is next taken up 
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in Section [5} This discussion begins with a review, in our Schwarzschild-form variables, of 
the classic analyses of delta-function sources of Refs mm- In the full nonlinear higher- 
derivative theory, exact solutions are not known and so one must use perturbative ex¬ 
pansions within the various Frobenius indicial families in order to analyse coupling to a 
delta-function shell source. Identifying the ‘vacuum’ with the non-singular (0, 0) indicial 
family and requiring this to be the solution type occurring inside a shell source, we find 
that only the (2, 2) indicial family has the correct number of free parameters required to 
match the various continuity and jump conditions needed across the delta-function source. 

Linking what happens near the origin to the behaviour of solutions near spatial infinity 
becomes the next issue to be considered. In Section [6j we generalise the result of Ref. 
m to show that the Ricci scalar in a portion of spacetime with Minkowski signature must 
vanish for any asymptotically-flat solution in the (0, 0) or (1,-1) indicial families. For these 
indicial families, this result obtains regardless of whether one considers a solution with a 
horizon at some intermediate radius ro, as in Ref. [14] . or considers a solution without a 
horizon. Requiring R = 0 correspondingly reduces the number of free parameters by one in 
each of these (0, 0) or (1, —1) indicial cases. 

For the traceless part of the higher-derivative field equation, the situation is complicated 
by errors made in the analysis of Ref. [H], as reported previously in Ref. m- Details of 
the corrected calculation are given here in Appendix [Cj One consequently does not have a 
straightforward way to prove a complete no-hair theorem setting the full Ricci tensor to zero 
in the (0, 0) or (1, —1) cases. However, for asymptotically-flat solutions with a horizon, one 
can still use linearised perturbation theory starting from the Schwarzschild solution. First, 
in Section [7J we use a Frobenius analysis about the horizon to show that such solutions, 
subject also to the requirement of a vanishing Ricci scalar as found in Section [6j have 
just three free parameters. This parameter count identifies the corresponding solution 
family with the indicial (1,-1) family near the origin, subject also to the requirement 
of asymptotic flatness as r —> oo and hence requiring also a vanishing Ricci scalar. The 
classic Schwarzschild solution is of course itself a member of this family, with just two free 
parameters (corresponding to the horizon radius and to a trivial time-rescaling parameter). 
Accordingly, the higher-derivative theory admits just one ‘non-Schwarzschild’ parameter 
controlling deviations from the Schwarzschild solution. 

Deriving a perturbative no-hair theorem for solutions expanded to linear order in the 
non-Schwarzschild parameter is then carried out in Section [HJ For a given horizon radius 
ro, and for asymptotically-flat solutions treated to linear order in the non-Schwarzschild 
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parameter, one finds that there is a range of small values of 01/(7r/j) for which the only 
asymptotically-flat solution with a horizon is the Schwarzschild solution itself. The precise 
range of such values depends on optimisation details of the linearised no-hair theorem, 
but the range boundary turns out to be quite near the phase bifurcation point for non- 
Schwarzschild black hole solutions found numerically in Ref. m- It may turn out that 
the range boundary for the linearised no-hair theorem and the black-hole phase bifurcation 
point actually coincide. 

In Section [9j we consider the more difficult question of what happens more generally in 
between the origin and spatial infinity. Owing to the complexity of the field equations, this 
can only be approached by numerical methods. One family of solutions at the origin that 
can mesh with the structures found at spatial infinity is the (2, 2) family. This agrees with 
numerical results found in Ref. m for the specific theory with m -2 = m 0 (i.e. the theory 
with a = 3/3). Generically, such (2,2) solutions have six free parameters at the origin, and 
six parameters at infinity, of which two combinations in each set must be adjusted in order 
to kill rising exponential behaviour from the spin-two and spin-zero sectors of the theory, 
thus leaving a four free-parameter set at the origin corresponding to the four parameters 
occurring at asymptotically-flat spatial infinity. We give another illustration of such a 
solution for the 7 R — aC 2 theory (i.e. with (3 = 0), in which the equations simplify owing 
to the absence of the spin-zero mode, then displaying a restricted match between a three- 
parameter set at the origin and a three-parameter set at asymptotically-flat spatial infinity. 
These (2,2) solutions cannot have horizons, since we have established in Sections [6] and [ 7 ] 
that asymptotically-flat solutions with horizons must belong to the (1,-1) family. Instead, 
asymptotically-flat solutions displaying Yukawa massive corrections at spatial infinity track 
closely to the Schwarzschild solution far out from the radius where the Schwarzschild horizon 
would have been, but they then begin to differ strongly from Schwarzschild as one comes in 
toward smaller radii, failing to have a horizon but matching instead onto the (2, 2) indicial 
family of solutions near the origin, and displaying a naked singularity. 

Section [9] also considers the structure of the (1,-1) solution family. The conclusion 
one draws from the linearised no-hair theorem of Section [8] is that the Schwarzschild solu¬ 
tion is generally isolated within the family of (1,-1) asymptotically-free solutions with a 
horizon, except for values of ( = ck/( 7^) located above a certain value C max , which presum¬ 
ably may be identified with the black-hole phase bifurcation point. Below this bifurcation 
point, perturbation in the single non-Schwarzschild parameter away from the Schwarzschild 
solution within the (1,-1) family initially can only lead to non-asymptotically-flat solu- 
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tions. The linearised no-hair theorem is thus in full agreement with the conclusions found 
numerically in Ref. m- The (1,-1) family naturally contains the Schwarzschild solution 
itself, with non-asymptotically-flat solutions generally occurring nearby as one adjusts the 
non-Schwarzschild parameter. As found in Ref. [T5] . however, there do exist additional 
asymptotically-flat (1,-1) solutions with a horizon that in general must be distinctly sep¬ 
arated from the Schwarzschild solution in the value of the non-Schwarzschild parameter. 
Such non-Schwarzschild solutions occur in the £ < Cmax range for which the perturbative 
no-hair theorem is applicable. This other branch of asymptotically-flat black-hole solutions 
accordingly exists for horizon radii ro greater than a certain value r™ 111 . At spatial infin¬ 
ity, such non-Schwarzschild black holes have a f e ~ m2r Yukawa correction to the goo metric 
component in addition to the 2 M/r Newtonian term, where M > 0 is the ADM mass. 

Numerical study using the shooting method for the horizonless (2, 2) solutions and for 
the non-Schwarzschild black-hole solutions reveals another feature of the overall spherically- 
synnnetric and asymptotically-flat solution space. Such solutions, with well-understood 
behaviours in each (small r and large r) asymptotic region of the radial coordinate r, 
appear to lie on separatrices between numerically found solutions with differing kinds of 
divergent behaviour. The implications of this separatrix structure for the overall solution 
space remain to be more fully understood. 

Another type of asymptotically-flat solution that emerges from numerical study may 
be described as a ‘wormhole’. In such a solution, which we also discuss in Section [9] for 
the (5 = 0 theory, the inverse of the A = g rr component of the metric goes to zero but 
the — B = gu component does not. General Z2 symmetric solutions of this type are highly 
constrained, with only two free parameters: the trivial time-rescaling parameter and the 
radius ro at which 1 /A vanishes. Numerical results show that such solutions can achieve 
asymptotic flatness at spatial infinity only for a particular value of ro, which is presumably 
related to the y / 2a/7 length scale. The Z2 symmetric wormhole solution is also found to 
lie on a separatrix lying between less regular solutions. 


In the Conclusion (Section 10) we give a brief discussion of some possible physical 
implications of our results. Clearly, the physical relevance of the present analysis depends 
upon fully accepting the implications of the higher-derivative terms in the field equations 
for the theory’s solutions, instead of simply considering their effects as perturbations on 
the second-order Einstein theory. It is equally important that there be at least some range 
of energy/length scales at which the fourth-order terms dominate, without their being 
swamped by the effects of yet higher-order terms. Given such assumptions, we comment 
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on stability questions for the various spherically-symmetric solutions, and on the possibility 
of phases in which the classic Schwarzschild solution might itself turn out to be the most 
stable. 


2 Fourth-order equations of motion 


The equations of motion derived from the action (1.1) are 
1 51 


H[lv ' V=gSg^ 


( 2 . 1 a) 


= 7 ^ R^ - ^g^uRj + ^ (a - 3i/3) Vfj.VuR - 2 aUR^ + * (a + 6(3 ) g^UR 

—AaR^R^vx + 2^/3 + RR^w + ^2aR' lX R 71 \ — ^/3 + -oA ^ 


1. 


— -T 

— 2 ’ 

satisfying a generalised Bianchi identity 


V" H ^ = 0 


and with trace 


= 6/3Di? - 1 R = -TJ 1 , 


( 2 . 1 b) 


( 2 . 2 ) 


(2.3) 


M ' I 2 ^ 

which is of fourth-order in derivatives of the metric for /3 7 ^ 0 and of second-order for (3 = 6. 
In fact the (3 = 6 (Einstein-Weyl) theory will turn out to be of particular interest to us. Note 
that in the Einstein-Weyl theory we can identify the more desirable sign of a by linearising 
around a Minkowski background. Writing g )w = + h^, we find 


- |aD(D - -) h^ = 0 


(2.4) 


And so a. > 0 is required for the absence of tachyonic instabilities. 


From the study of the linearised limit of (2.1) about flat spacetime in [3], one knows 
that there are massive spin-two and spin-zero excitations with masses 


m o : = 


m n : = 


2a ’ 
J_ 
6/3 ’ 


so one notes that 


H£ = 6/3 (□ — m 0 2 ) R ■ 


(2.5a) 

(2.5b) 

( 2 . 6 ) 
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When considering spherically-symmetric static solutions, we may take the metric to 
have the Schwarzschild form 


ds 2 = —B(r) dt 2 + A{r) dr 2 + r 2 dd 2 + r 2 sin 2 6 dcj) 2 , 


(2.7) 


and we shall look for source-free solutions with = 0 in the bulk of spacetimej^] With 
this ansatz there are only two independent equations of motion, for A{r) and B(r). The H 
field equation tensor takes the form 


H IW = 


(H tt (r) 0 0 

0 H rr (?’) 0 

0 0 H ee (r) 

\ 0 0 0 


\ 


( 2 . 8 ) 


Hoo(r) sin 2 0/ 


the components of which are related (for 6 = 5) by the r component of the Bianchi identity 


( 2 . 2 ): 


H r r 

~A 


2 H rr B'H rr 2 Hqq B'H tt _ 

+ Ar + 2AB + 2 B 2 = ' 


(2.9) 


Accordingly, when provided with a stress-tensor source, the system is described by just 
two independent equations 


Hu = -T, 


tt : 


Hrr — i r 


(2.10a) 

(2.10b) 


If the metric (1.3) is substituted into the Lagrangian before performing the variation, the 


resulting equations of motion are found to be equivalent to the set (2.10), i.e. the truncation 
to the static spherically-symmetric case is a consistent truncation: 

si _ V=g 


<5A 

51 _ 

6B 


A 2 H rr , 


V^9 

B 2 


Hr 


(2.11a) 

(2.11b) 


Consistency of truncation to the metric form (1.3) is guaranteed in the usual fashion because 


one is truncating to the invariant sector under a group action - in this case spatial rotations 
[IB]. It should be emphasised at this point that in this paper we are not making any 
additional simplifying truncations such as setting AB = constant. Imposing such additional 


5 We shall address the issue of delta-function sources for solutions in Sections |4j and JHJ 
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conditions certainly makes solution of the equations greatly simpler, but it also severely 
restricts the corresponding solution set j^] 

From here onwards, unless otherwise stated, we will be solving the source-free equations 
for r > 0 

H tt = 0, (2.12a) 

H rr = 0. (2.12b) 


2.1 Differential order 

2 . 1.1 /3/0 


To find the differential order of these coupled equations, note that Hu is a function of 
A) 3 )(r), B^ 3 \r), B^\r) and lower-order derivatives^ and H rr is a function of A"(r ), 5"(r), B^ 3 \r) 
and lower-order derivatives. Let us now analyse the differential order of these equations. 

Note that a = 0 and (3 = 0 are special cases of different differential order, so we shall first 
look at the generic case a / 0, a / 3/3, (3 / 0. We define 


A 2 (r(a - 3/3)5' - 2(a + 6/3)B) 2 
x ((a - 3/3)5 (2 rBA' (r(a - 3/3)5' - 2(a + 6/3)5) 

+A (—r 2 (a - 3(3)B' 2 - 4 r(a - 3/3)55' + 12(a + 6/3)5 2 ) ) , 

Y( r ) = 2r(a — 3(3)B 2 

K J A (2(a + 6/3)5 - r(a - 3/3)5') ' 

It is clear that the equations of motion are equivalent to the pair 


(2.13a) 

(2.13b) 


0 = H rr , (2.14a) 

0 = H a - X{r)H rr -Y{r)d r H rr , (2.14b) 


the first of which is of third order in B and of second order in A, and the second of which 
is of third order in A and of second order in B. Details are given in Appendix [Aj This 
reduction in order leads us to expect a total of six free parameters in the solution. This can 
be more clearly seen by eliminating 5(r) to get an equation of sixth order in A(r) alone; 
the detailed procedure is sketched in Appendix |A| 

6 Indeed, in m, among other cases, the pure C^, vpa C ,1 '' pa theory was considered subject to such a restricted 
ansatz, with the result that the only solution without conical singularities is just the classic Schwarzschild 
solution. The same simplifying restriction has been made recently in the analysis of spherically symmetric 
solutions for the pure R 2 theory in Ref. m,■ Our aim in the present paper is to explore the full set of 
spherically symmetric solutions without such a restriction. 

7 We denote derivatives of order > 3 by superscripts such as A ^. 
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2.1.2 (3 = 0 


In this special case, the massive scalar is absent. In the trace of the sourceless equations of 
motion (2.3) with = 0, two derivatives disappear and the equation simply states that 
R = 0. As this suggests, the total differential order in this case is reduced by two with 
respect to the (3 / 0 case. The equations of motion are then equivalent to the pair 


0 = H M 

u 11 fi i 

_ H rr 3rBA' -2A(rB' + B) + 2A 2 B 
°~ + 3 7 r 2 AB 



d r (Hp 


2 B — rB' 
O^rB 


(2.15a) 

(2.15b) 


the first of which is of second order in B and of first order in A, and the second of which is 
of second order in A and of first order in B. 

These two second-order equations imply that there are four free parameters in the 
solution for the (3 = 0 case. 


3 Solutions near the origin and Frobenius analysis 


Previously in Ref. [3], the asymptotic behaviour of solutions to the equations of motion 
was analysed near the origin, working to leading orders in r. Here we will solve expansions 
to several higher orders in r in order to improve our understanding of the parametric 
dependences of solutions. 

The two undetermined functions in the metric are expanded in Frobenius series in r as 


A{r) = a s r s + a s+ ir s+1 + a s+2 r s+2 + ... , 
B(r) = bt {r t + b t+1 r t+1 + b t+2 r t+2 + ...) , 


(3.1) 


where a s ,bt 7^ 0 are nonvanishing coefficients and s,t are indices yet to be determined. 


Substituting the series (3.1) into the equations of motion (2.12) and analysing the con¬ 
sistent possibilities for the (s, t) indices leads exclusively to three solution families |i.3J for 
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generic^] a / 0, /3 / 0: 


M) = (o, o) 

(s,t) = (1,-1) 

M = ( 2 , 2 ) 


(3.2) 


In each of these families, the equations can be solved order-by-order for the coefficients 
a n , b n . Some coefficients will be left undetermined in this process, corresponding to the free 
parameters of the system in each solution family. There will always be one free parameter 
in B{r ) corresponding to a trivial scaling of the time coordinate. 

Of course, in performing asymptotic analysis of this sort, an assumption is being made 


that Frobenius type expansions such as those of Eqs. (3.1) with integral steps in powers of 
r following the leading (r s , r l ) terms is adequate to capture all possible types of asymptotic 


behaviour for solutions to nonlinear equations such as (2.12). One might worry about the 


inclusion of terms such as exp (c/r p ) times a Frobenius series, or of terms involving powers 
of logarithms. In the case of linear systems of differential equations, one can deal with 
such possibilities on the basis of general theorems about equation systems with regular or 
irregular singular points of various ranks, but a suspicion could remain that this might not 
capture the full complexity of solutions to systems such as ( 2.12| ) (written out in full detail 
in Appendix [A]). All we can say to dispel such concerns is that we have explicitly tried many 
such exotic possibilities and the only consistent leading asymptotic behaviours that we have 


found are those shown in (3.2). 


3.1 Free-parameter counts in each of the near-origin solution families 


We have expanded and solved the equations of motion (2.12) to at least twelve orders in r. 


In each family, all free parameters have appeared by the fourth order at the latest, and after 
that each new order brings two new parameters and two new constraints. The resulting 
free-parameter counts are given in Table [lj 


We neglect here families of solutions appearing only at special values of a > 3/3 > 0: 


t - 2 


= s € 


(s 2 + 2 s + 2) 2 
a = i - t- —3/3 


and two families of solutions for a = 0: 


4 + 2t + r 
4 + t 


= se 


or (s, t) = (0,1) (with 1 free parameter). 
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Table 1: Free parameter counts for the indicial solution families. 


(s,t) solution family 

number of free parameters 

a choice of free parameters 

( 0 , 0 ) 

3 

bo, CL2, b 2 

( 1 ,- 1 ) 

4 

ai, b-i, a 4 , b 2 

( 2 , 2 ) 

6 

02 , b 2 , 63 , 64 , 05 , 65 


3.1.1 The (0,0) family 

The first few terms in the (0, 0) family are: 


A(r) 


B(r) 

b 0 


1 + o 2 r 2 

4 a 2 6 o (&o7(2a + 3/3) — 36a/36 2 ) + 18a|/36Q(10a + 3/3) — 2 b 2 ( 607 ( 0 : — 3/3) + 9/36 2 (2a + 3/3)) 
r 180a/36g 

+0(r 6 ), (3.3a) 

1 + b 2 r 2 

r 4 (54a|/3 2 + a 2 (-07 + 108a/36 2 + 3/3y) + b 2 ( 7(0 + 6/3) + 54/36 2 (2a - /?))) 

360a/3 

+0(r 6 ). (3.3b) 


The three-parameter (0, 0) solution is the natural ‘vacuum’ solution family of the higher- 
derivative theory, comparable to the two-parameter spatially homogeneous flat space solu¬ 
tion in Einstein theory. The Riemann curvature tensor R a bcd referred to an orthonormal 
frame is nonsingular as r —> 0 for this solution. 


3.1.2 The (1,-1) family 

The first few terms in the (1,-1) family are: 


A(r) = air — a\r 2 + a 3 r 3 + a 4 r 4 — —( ai (^ a i b 2 + 19 af + 3504 )) 

+-^a 2 r 6 (21ai6 2 + 101a 4 + 141a 4 ) + 0(r 7 ), (3.4a) 

= 1 + ai + b 2 r 2 + -^r 3 (ai& 2 + a 4 + a 4 ) 
o_i r 16 

-^3r 4 (ai (ai& 2 + af + a 4 )) + 0(r 5 ). (3.4b) 


The (1,-1) family is clearly the family that contains the classic Schwarzschild solution of 
Einstein theory. The Schwarzschild solution is obviously a solution of the higher-derivative 


theory’s field equations because every term in (2.1) contains R or R. At the origin, the 
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(1,-1) indicial structure gives rise to a curvature singularity, with R /llip(T R IJ ' l ' p(T going like 
r -6 as r —> 0 0 . 


3.1.3 The (2, 2) family 

The first few terms in the (2, 2) family are: 


A(r) = a 2 r 2 + a 2 b 3 r 3 - (2a 2 + 6 2 - 864 ) + a 5 r 5 


+ 


( 


- 12 a 2 a\ - 2a\ ( 6 § (a 2 - 603a/3 - 252/3 2 ) + 27a (20/36 4 + 7 )) 


1296a/3 \ 

+a 2 (-16a 2 + 1413a/3 - 72/3 2 ) + 2b A b\ (19a 2 - 2223a/3 + 180/5 2 ) 


62 


—366563 (a 2 + 45/3 2 ) + 12a6|(a + 162/3)J + 324 as/ 363 ( 7 a + 3/3) 
+ 0 (r 7 ), 

= r 2 + & 3 r 3 + 6 4 r 4 + 6 sr 5 
„6 


+ 


216aa 2 


— 12aa 3 + a 2 (l46 2 (2a + 3/3) — 24a6 4 ) 


(3.5a) 


+a 2 (26^(67a - 3/3) + 26 4 6f (15/3 - 227a) + 456 5 6 3 (7a - 3/3) + 180a6^) + 27a 5 6 3 (a + 3/3) J 

+0(r 7 ). (3.5b) 

The (2, 2) indicial structure at the origin gives rise to a curvature singularity for this 
solution family, with Rpv pa R p ' l ' pa going like r ~ 8 as r —> 0 [3j. 


4 Coupling to sources in the linearised theory 


For asymptotically-flat solutions, the weakening fields as r —> 00 can reliably be analysed 


using the linearised limit of the field equations (2.1). We now set the stage for our later 


discussion of source coupling in the full nonlinear theory by studying coupling to sources in 
the linearised theory, expanding somewhat the discussion given in Ref. [3j. We first need to 
identify the vacuum solutions that can occur inside a shell source. 


4.1 Source-free solutions 

In |3], the linearised equations were solved for 

A = 1 + W(r) + 0(W 2 ), 

= 1 + V(r) + 0(V 2 ). 


B 
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(4.1a) 

(4.1b) 







Solving the linearised source-free equations (2.12) for r > 0 in this limit yields the general 
solution 


V = 
W = 


n , c 2 ,0 , Co-e~ m ° r C 0+ e mor C 2 -e~ m * r C 2+ e m * r 

O -| | | | | 


r 

O ~mor 


O mor 


f' p U(>|jl /jll 

—— + Co- -(1 + m 0 r) + C 0 + -(1 - m 0 r) 


-~C 2 -- 

2 


0 —rri2r 


1 


0 rri2r 


-(1 + m 2 r ) - -C 2+ - 

r 2 r 


(1 — m 2 r) . 


(4.2a) 


(4.2b) 


Note here that for a < 0 or f3 < 0 one has pure imaginary m 2 or pure imaginary mo, respec¬ 
tively. For pure imaginary masses m, = i /i, the solution has, instead of real exponentials, 
suppressed oscillating terms like C lyS £ sin(^r) and C 1yC ^ cos (mr) in V and W. However, 
W also has non-suppressed oscillations behaving like Ci >s sinQujr) and Ci fC cos(/rjr). This 
precludes asymptotic flatness at spatial infinity unless both of the constants Ci vanish, i.e. 
unless the corresponding metric solution is strictly flat. Accordingly, we limit our consider¬ 
ation to cases where m 2 > 0 and mg > 0. 


The linearised solution (4.2) clearly shows the existence in general of six free parameters, 


noting that the free parameter C corresponds to a trivial rescaling of the time coordinate. 


As one approaches the origin in the linearised solution (4.2) for generic values of the free 
parameters, the Cartesian-coordinate form of the linearised curvature tensor R]^ cd has lead- 


3 


mg r ” singular terms in the r —> 0 limit. The linearised Ricci scalar for the solution (4.2) 
is 


R = - 


3C 0 _mge- mor 3 C 0+ mle mor 


r 


(4.3) 


and so has leading r -1 behaviour for generic values of the free parameters. 

Note that a = 0 or (3 = 0 or a = 3/3 are special cases in which m 2 or mo vanish or 
coincide. In the following, we shall proceed for the generic case a / 0, /3 / 0, a / 3/3. 

4.2 True linearised vacuum 


When the general linearised solution (4.2) is extended all the way to the origin at r = 0, H, 


flV 


must in general involve S 3 ^) sources. The true vacuum solution without such delta-function 
sources is the restricted three-parameter solution family satisfying the vacuum constraints 


C 2 .o — C 2 - + C 2+ — Co- + Cq+ — 0 . 


(4.4) 
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Subject to these constraints, one finds the true linearised vacuum family 
2 sinh (m 0 r) 2 sinh(m 2 r) 

^vac — ^ I ^0+ i O2+ 

r r 

/ sinh(m 0 r) A „ / sinh(m 2 r) 

IVyac = 2C 0 +-m 0 cosh(m 0 r) - C 2+ -m 2 cosh(m 2 r) 


(4.5a) 


(4.5b) 


A consequence of the vacuum constraints (4.4) is that the metric (4.1) for the true linearised 


vacuum family (4.5) is nonsingular at the r = 0 origin. This nonsingularity extends as well 
to all components of the linearised curvature tensor i?|T. and in particular one can see 


from (4.3) using (4.4) that the linearised Ricci scalar is nonsingular at the origin for the 


true vacuum solution (4.5). 


It is appropriate to distinguish the true vacuum solution (4.5), with its nonsingular 


curvature, from other members of the general solution family (4.2) that happen to have a 
nonsingular metric as r —> 0. Unlike the situation in linearised Einstein theory, where the 
only spherically symmetric solution with a nonsingular metric at r = 0 is simply flat space, 
with correspondingly vanishing curvature, in the linearised version of the higher-derivative 


theory ( 1 . 1 ) the family of nonsingular-metric solutions turns out to be wider than just the 


vacuum solution (4.5). This wider class of nonsingular-metric solutions includes also the 


solution for a point delta-function source, which we consider next. 


4.3 Source examples for the linearised theory 
4.3.1 Point source 

In [3j the stress-tensor of a static point mass at the origin was considered: 

T„„ = 5y v M5\x ). (4.6) 

With this source, the solution to the linearised equations of motion is the vacuum solution 
plus an asymptotically-flat matter part: 

M 


V(r) = C- 


0 


—mor 


— 4e 


—m^T 


+ 3) , 


247T7r 

W(r) = -((l + m 0 r)e- mor + 2 (l + m 2 r)e- m 2 r -3) 

247T7r 


(4.7a) 

(4.7b) 


indicating that one needs 7 = ^ in order to agree with the Schwarzschild solution 

in the limit where m,Q and m 2 tend to infinity. 


As one can see from (4.7), V and W are actually nonsingular as r —> 0. As is clear 


from the need for the source (4.6), however, (4.7) cannot be considered a true vacuum 
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solution. This conclusion is reinforced by consideration of the curvature R^ c d as r ~^ 0; 
whose components have leading r~ l singularity, and for which the Ricci scalar is given by 


R = 


M 

8777 r 


777-0 e 


—rriQr 


(4.8) 


which has r _1 behaviour as r —> 0. Note also that in the weak-field regime with 7770 finite 
{i.e. for /l / 0) a solution with a point source at r = 0 always has R 7 ^ 0 at any nonzero 
value of r. 


4.3.2 Shell source 


To illustrate the effects of extended sources in this theory with no Birkhoff theorem, let 
us now solve again in the linearised theory for the fields produced by various sources of 
nonzero size. 

First take as source a thin spherical shell of radius t 

M 


T u = 


4lX 


S(r-£), 


T rr = 0 . 


(4.9a) 

(4.9b) 


From the linearised = 0 condition, we have 


T e0 = O + O((W,R) 2 ). 


(4.10) 


For r < t we use the vacuum solution (4.5): 


Vin = D - 


2T>o- sinh^TroT’) 2D2- sinh(r7727’) 


Win = —2T>0- ( Sin * 1 ( m ° T ’) _ mo C osh(777 0 7‘) ) + D 2 - 


sinh(r 77 - 27 ’) 


(4.11a) 

— 7772 COsh(7?727’) 

(4.11b) 


and for r > l we use the source-free solution (4.2a) for T4 u t and (4.2b) for W out with the 


rising exponentials suppressed in order to achieve asymptotic flatness. 

For a/0 and /3 7 ^ 0, one finds that V'"{r) and W"(r) are discontinuous at the location 
of the shell, and the solution is 

T M (1 1\ e ~ m 2 r M sinh( 7772 ^) e _m ° r M sinh(? 77 oT) ,, iri . 

Rout = £> + — (-?--)+ - 7 --^- ^ir 1 , (4.12a) 


8777 \£ r 


r 6777 7772 ^ 


r 247T7 TTIqI 




247T7r 


mod 


7772-^ 


(4.12b) 
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for which the Ricci scalar is 


Rout, - 


M 


2 —mor 


m 0 e 


sinh(mo^) 


(4.13) 


87T7 r u m^l 

Thus, just as for the point source, we find in the weak-field regime with mo finite ( i.e. for 
[3 7 ^ 0) that the solution with a shell source always has R ^ 0 at any nonzero value of r. 

Note that in the limit l —> 0 of the shell-source solution, the expressions for V out , W ou t 
and R correctly tend to those of the point source. 


4.3.3 Balloon source 

Now let us take as source a stress tensor with internal pressure, again expanding upon 
results given in [3]: 

0 0 0 \ 

0 PQ(£-r) 0 0 

0 0 Tqq 0 

y 0 0 0 Tqq sin 2 9 

where 0(r) is a Heaviside theta function^] In order to satisfy the linearised conservation 
condition for T /w , we need to have 


T — 

fits — 


(4.14) 


T ee = Pr 2 Q(l - r) - ^Pr 3 S(l - r ). 

Solving the system with this source, we find the interior r < l solution 

. 2T>o-sinh (mor) 2L>2- sinh (mor) r 2 (4irl 3 P + M) 

Vi n (r) =--- v ; ----— PP + D + v y 


VTi n (r) = Do- ( 2mo cosh (mor) — 
l sinh (m 2 r) 


r 

2 sinh 


(- m 0 r) \ 


167T7 P 


+D 2 - 


— m 2 cosh (m 2 r) I + 


Mr 2 


(4.15) 


(4.16a) 


(4.16b) 


r j 87T'yl 3 ’ 

where the three vacuum constraints 0 = D 2 ^ = D 2 - + D 2+ = Do- + Dq+ have again been 
used to ensure a pure-vacuum r < £ internal solution without source. (Note that the M, P 
source terms are proportional to r 2 and do not affect these constraint requirements.) At 
r = l there are 5 continuity conditions, for V, V ’, V". W, W', and two step conditions 

lP{a + 6/3) 


Ct (i+)=K'(£-) + 


WUU)=W^(£-)- 


36 a/3 

l 2 P(a - 3(3) 
36 a/3 


(4.17a) 

(4.17b) 


9 « 


0 (r) is capitalised in order not to be confused with the angular coordinate 8. 
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Note, however, that of these seven continuity and step conditions, only six are independent 


( c.f. the six free parameters expected from the differential order analysis in Section 2.1). 
This is a general situation, and will be important for us when we consider such source 
couplings in the full nonlinear theory. 

Implementing the continuity and step conditions, we obtain the asymptotically flat 
solution with a balloon source: 

„2 


Mn =D + 


1 

487T7 1 


3M P £ + 2 [3 (1 + mol)M 0 - 4v d 3 P] smh ( m ° r ) e ~™oi 
l z L J mor 


- 8 [3 (1 + m 2 l) Mi + 2vr/ 3 P] Slllll ( m2r ) e ~m 2 i 


W in = 2 ^ ( 3M^ + [3 (1 + mol) M 0 - 4tt/ 3 P] 


m 2 r 
sinh(mor) 


(4.18a) 


mor 


— cosh(rao'r) 


0 —mol 


sinh(m2r) 


+ 2 [3 (1 + m 2 l) M 2 + 2nl 3 P] 

[ m 2 r 

bout = D + ——- (2Mq — 8M 2 + 3 M + 4tt1 3 P) — — - 

lo7ry£ v ' 8n r yr 


— cosh(m2r) 


D —m,2l 


+ 


,-m 0 r 


24:7r-yr 

f,-m 2 r 


Q'K'yr 


3 Mo 


3M 2 


sinh(moO 

mol 

si nil (to 2 0 
m 2 l 


— cosh (mo l) 

— cosh(m 2 l) 


_ 4 ^ 3p sinh(TOoO \ 
mol ) 


+ 2t xl 6 P 


3 sinh (m 2 l)\ 


Wo„ t =- 


M 


8'JT'yr 

e -m 0 r _|_ mor ) 

247 T 7 r 

e -m 2 r( i _|_ m2r ^ 


127T7r 


3M 0 

3M 2 


sinh(moO 

mol 

sinh(m2Z) 

m 2 l 


— cosIi(tooZ) 

— cosh(m2/) 


+ 2t tPP 


3 sinh(TO2/) \ 


m 2 l 


(4.18b) 


m 2 l J 


4 ^ 3p sinh(TOoO\ 

m 0 l J 


(4.18c) 


J (4.18d) 


where we have used the following notation for source-parameter combinations 


M P 

Mo 


M + 4t d 3 P 
M - 4t t1 3 P 


Z 2 mg 


(4.19) 


, r M + 2nl 3 P 

M2 := Pm\ ■ 

The main point to take away from this analysis of the linearised solutions is that the 
general six-parameter solution, constrained by two requirements of vanishing coefficients 
for the rising-exponential Yukawa terms as r —> oo, has a remaining essential dependence 
on four parameters. One of these is adjustable by rescaling of the time coordinate t (cor¬ 
responding to the additive parameter D above), and will be fixed by the requirement of 
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having an asymptotic Minkowski metric as r —> oo. The other three parameters will be fixed 
by details of the source, as displayed in the balloon-source solution by the dependence on 
£, M and P. This multi-parameter dependence clearly illustrates the absence of a Birkhoff 
theorem for the higher-derivative gravity theory. One needs to start with the full six- 
parameter generic solution in order to arrange a successful coupling of the higher-derivative 
theory to a standard matter source, exemplified here by these various delta-function source 
constructions. 


5 Shell sources in the full nonlinear theory 


We now progress to studying matter coupling in the full nonlinear theory. Unlike the 
situation in general relativity, where the Schwarzschild solution is known in closed form, we 
have no such luxury in the higher-derivative gravity theory. So we need to be careful in 
handling the continuity and step matching conditions for solutions known only from series 
expansions such as those given in Section |3j What we wish to establish is which of the 
three families (s,t) = (0,0), (1, —1) or (2,2) can couple acceptably to an ordinary matter 
stress-tensor source. The key to this will be the parameter counts that we found in Section 

El 


For simple models of matter coupling, we again consider distributional sources. As has 
long been clear m in general relativity, however, the only sensible delta-function sources in 
generally covariant theories are sources of spatial codimension one. So we do not consider a 
point source as in the linearised theory. Instead, the simplest source that we can consider in 
the full nonlinear theory is a thin spherical shell of radius £, which can be compared to the 


discussion given for the linearised theory in Section 4.3.2 This shell source has a conserved 
stress tensor 


T — 

/IV - 


(T tt 
0 
0 
\0 


0 

0 

Tee 


\ 


(5.1) 


0 Tqq sin 2 (0)y 


where, as in (4.9), 


T tt = 


M 
4-71T 2 


S(r-£), 


T rr = 0. 


(5.2a) 

(5.2b) 
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The condition V ll T lw = (0, V^T^, 0,0) = 0 requires 


Tee 

The equations of motion (2.1) expand 

_ r*B'T tt 

4 B 2 ' 

schematically as 

(5.3) 

H tt = ~ B^+ 

~ A( 3 )+ ~ H (3 ) + ... , 

(5.4a) 

H rr = ~ H (3) + 

~ A"+ ~ B" + ... , 

(5.4b) 

Bee = ~ B^p 

~ A (3) + ~ H (3) + ... , 

(5.4c) 

suggesting that we should consider 

B ( 4 ) 

~ 5 + 0 , 

(5.5a) 

A (S) 

~ 5 + 0 , 

(5.5b) 

Hi 3 ) 

~ 0, 

(5.5c) 

A" 

~ 0. 

(5.5d) 


Then A, A', B , B ', B" will be continuous at r = l. while A" has a step of size 


A" (c \ _ A" (f \ - M ~ 3/3 ^- g/ ~ 2 ( a + 

out! +) ini -) 87r£ A 36a/3 


(5.6) 


We leave to Appendix [B] a detailed discussion of how to arrange a satisfactory series so¬ 
lution of these matching conditions in the higher-derivative theory. However, the important 
part of the result is easily seen by a parameter-counting argument as follows. The region 
interior to the shell is described by the vacuum solution of the nonlinear theory, which is 


the (0,0) family as discussed in Section 3.1.1 The (0,0) vacuum solutions of the nonlinear 
theory admit three free parameters, as shown in Table [l] given in Section [3} The coupling 
to the source constitutes six continuity and matching conditions. The region exterior to the 
shell is also source-free, so will be described by one of the source-free solutions of Section [3j 
Finally there are also two conditions at infinity that need to be imposed in order to ensure 
asymptotic flatness, analogous to the elimination of the rising-exponential Yukawa terms in 
the linearised theory. After applying all these constraints we expect one final free parameter 
to be the adjustable parameter in the B function, corresponding to the asymptotic value of 
goo; which needs to equal —1 in order to have asymptotic Minkowski space as r —)• oo. The 
remaining structure of the solution will be determined by the details of the source, deter¬ 
mined by the two parameters M and i in the case of the simple shell delta-function source, 
and determined by the three parameters M, £ and P in the case of a balloon-type source 
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as in Section [4.3.3 This is a total of 6+2+l=9 conditions for asymptotically Minkowski 
shell-coupled solutions, with 3 free parameters from the (0,0) solution inside the shell, so 
the solution outside the shell must have 6 (or more) free parameters in order to be able to 
satisfy the constraints. Thus the exterior solution to a matter shell must be of the (2, 2) 
family, which has precisely 6 free parameters. This is similar to the linearised source cou¬ 
plings where the exterior solutions had all of their 6 free parameters fixed by the parameters 
of the source, i.e. the generic 6-parameter solution can be placed outside the source and 
can satisfy the necessary constraints, but a constrained exterior solution could not. Note 
that we assume that the 9 continuity, matching and asymptotic flatness conditions are all 
independent. In the linearised theory, one can verify that the conditions are indeed inde¬ 
pendent, but strictly speaking only a closed-form solution could confirm that the same is 
true in the full nonlinear theory. 

On the other hand, trying to arrange the coupling of a delta-function shell source to 
an exterior (0,0) or (1,-1) family solution will not work, for the simple reason that their 
numbers of exterior free parameters (three or four, respectively) are not sufficient to satisfy 
all the nine continuity, step, asymptotic flatness and asymptotic Minkowskian requirements. 
Thus, an asymptotically-flat and asymptotically-Minkowskian solution coupled to a shell 
delta-function source can only be of the (2, 2) family. 

This contrasts with general relativity, where the (0, 0) vacuum family has 1 free param¬ 
eter, the (1,-1) family has two free parameters, and the (2,2) family does not exist. The 
shell coupling gives two conditions, gmij" ~> oo) = —1 gives one condition, and there are 
no conditions needed to remove asymptotically non-flat terms. The solution exterior to a 
shell source is the Schwarzschild (1,-1) family, with structure determined by M and the 
asymptotic Minkowskian condition, but is independent of i. 


6 Trace-equation no-hair theorem 

Having established that the solution that couples correctly to an ordinary stress-tensor 
source is of the (2,2) family, we now proceed to investigate the consequences of the field 
equations in the (0,0), (1,-1) and (2,2) solution families without regard to sources. We 
will be particularly interested in the consequences of boundary conditions at a putative 
horizon or at spatial infinity. Useful tools to this end are a set of Lichnerowicz-type ‘no¬ 
hair’ theorems forcing the solution to share properties with the standard Schwarzschild 
solution under certain conditions. This topic was broached in Ref. [T3]. As noted in Ref. 
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[I5|. we agree in part with conclusions of that reference, namely the trace part of the no¬ 
hair theorem, as will be discussed in the following. Unfortunately, we do not agree with the 
contentions of El regarding the traceless part the higher-derivative equations of motion, 
which would have significantly simplified the analysis. We present in Appendix 0 our 
analysis of the general no-hair theorem, including an extension to include a cosmological 
constant. Despite our disagreement with Ref. |14j , we still can obtain important constraints 


on the solution families of the higher-derivative theory (1.1) using no-hair type arguments. 
In this section, we review the trace no-hair theorem of 


. In Section [7| a careful analysis 
of the parametric structure of solutions containing a horizon will be given, and in Section [8] 
these elements will be put together with a no-hair theorem for linearised deviations from the 
classic Schwarzschild solution. In Section [9j solutions with a horizon that discretely differ 
from the Schwarzschild solution will be discussed. The only family of static spherically- 
symmetric and asymptotically-flat solutions that couples properly to ordinary stress-tensor 
sources, i.e. the (2,2) family, cannot have a horizon. 

We now proceed to review the trace-equation no-hair theorem. We do this in a different 
style from that of Ref. El in that we present the no-hair trace-equation argument for static 
solutions in terms of a timelike dimensional reduction from four to three dimensions!^! 

The static four-dimensional metric can be written in the form 


ds 2 = — A 2 dt 2 + h a h dx a dx b 


( 6 . 1 ) 


where the spatial metric h a b is positive definite for flat space, and therefore, assuming 
asymptotic flatness, it is positive definite everywhere between infinity and a horizon at 
finite r (should one exist). Both h a b and A are assumed to be functions only of the three 
spatial coordinates x a . Let be the covariant derivative for the 4-metric g^ v , and let D a 
be the covariant derivative for the 3-metric h a b- It follows that 


□-R := = D a D a R + J (D a A) (D a R) , 


so the trace of the source-free equations of motion (2.3) can be written as 


1 


0 = H' 1 = 6/3 ( D a D a R + - (D a A) (D a R) - rnJR 


( 6 . 2 ) 


(6.3) 


For the theory with (5 = 0, one has R = 0 directly from the trace of the field equations 


(2.3), while for (5^0 the bracketed quantity in (6.3) is required to vanish. Subject to 


10 Such timelike dimensional reductions have proven to be powerful tools in classifying black-hole solutions 
in a variety of supersymmetric and non-supersymmetric contexts [I8l . 
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certain assumptions, this will still imply that R = 0. To see this, multiply the bracketed 
quantity in (6.3) by XR and integrate over the whole spatial 3-section 


0 = j Vh d 3 x[X R ( D a D a R ) + A R (D a A) (D a R) - m 0 2 A R 2 ] 

Preparing for an integration by parts, we rewrite this as 

0 = [ Vh d 3 x [D a (A R D a R ) - A ( D a R ) (D a R) - m 0 2 A R 2 ] 
Js 


(6.4) 


(6.5) 


The first term turns into a 2-dimensional integral over two boundaries: one at spatial 
infinity, and the other at some finite radius. The contribution from the boundary at spatial 
infinity vanishes subject to the assumption of asymptotic flatness. A sufficient (though not 
necessary) condition for the inner boundary term to vanish is satisfied if that boundary is a 
horizon. If we take the inner integration boundary to be such a horizon, then the boundary 
integral will be proportional to A|horizon, which vanishes, by definition. Therefore, for a 


solution with a horizon, (6.5) reduces to a 3-dimensional spatial integral over the sum of 


two negative semi-definite terms j^] Requiring this to vanish therefore implies that R and 
D a R separately vanish throughout the integration region, thus implying that R = 0. 

In a region where R = 0, the equations of motion become 


0 = H, 


^\R=0 


= -2 a (nR^ + 2RPR vp - 2V p V p R p - ^g pv R pa R i 


pa 


+ 7 R 


11V j 


( 6 . 6 ) 


which notably no longer have any dependence on (3. In the (3 = 0 case without sources, one 


automatically has H p = qi? = 0, and therefore (6.6) obtains everywhere. 


We have accordingly shown using (6.5) that if the boundary terms vanish on the bound¬ 


aries of a given spatial region, then the field equations in that region reduce to the (3 = 0 
case (corresponding to the Lagrangian density C ~ 'yR — aC flupa C pL ' pa ). It will be com¬ 
putationally advantageous in such situations to use the two independent (3 = 0 equations 
of motion, rather than the two (3 ^ 0 equations of motion together with R = 0 as a third 
condition. 


6.1 Implications for the three near-origin solution families 

We have anticipated above a main conclusion that the asymptotically-flat solutions to the 
higher-derivative theory 0 with normal matter coupling do not have a horizon. For a 


spacetime without horizon it is natural to extend the integration region in Equation (6.5) 
down to near the origin, r —> 0, where h a {, remains positive but we do not yet know the 


1 Recall that we are requiring throughout /3 > 0, so > 0. 
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behaviour of the boundary term at the inner boundary. An analysis of the equations of 
motion at small r can, however, tell us how this boundary term behaves. 


We now look again at (6.5) but now with the integration region taken to have an inner 


boundary located at r —> 0. We can use this discussion together with the small r behaviour 


of the boundary term to study the various implications of (6.5) for the structure of the 


(0,0), (1, —1) and (2,2) solution families. Accordingly, we calculate the boundary term in 


(6.5) simply taken as the dominant part 


Rd r R 


(where the other components of the boundary term vanish in all cases as r —> 0). To 
be precise about the boundary-term contribution, recall that the boundary term actually 


appears in (6.5) in the form 


Vh 


D a I v B (boundary term)“ 


d 3 x 


(6.7) 


which for a boundary term with vanishing 6,cj) components is equal to 


d r 


\Th \Z~B (boundary term) T 


d 3 £ = J d r 


\[AB r 2 sin($) (boundary term) r drd6d(f> 

( 6 . 8 ) 


We now analyse the consequences of (6.5) in the three near-origin solution families. 


Family 

Rd r R 

\J~hB Rd r R 

(0,0) 

f (a 2 - b 2 ) 2 r + 0(r 3 ) 

~ 0(r 3 ) 

(1,-1) 

4 ° 4 ) r + 0(r 3 ) 

~ 0(r 3 ) 

(2,2) 

(a 2 (l4a 2 & 3 —2ft| + 10ft 4 b3—45& 5 )+27a 5 ) 2 5 n ( r -4) 

9a 5 2 r ) 

~ 0(r -1 ) 


For a spacetime with no horizon, we choose the integration region of (6.5) to extend from 
the origin to infinity. The boundary at infinity gives zero by the assumption of asymptotic 
flatness. For the (0,0) and (1,-1) families, as the inner integration boundary is taken 
towards the origin, r —> 0, the inner boundary terms also tend to zero, and we consequently 
learn that if there is no source between the origin and infinity then one must have R = 0 
throughout spacetime. The (0,0) family contains an R = 0 solution, as does the (1,-1) 
family. For the (2,2) family, the boundary term blows up as r —> 0 and one can make 
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no conclusion about the necessary vanishing of R. Note, however, that the (2,2) family 
also does contain an R = 0 solution, obtained by applying two constraints to its six free 
parameters. Compare this to the expression for the Ricci scalar in the linearised theory 


(4.3) and to the analysis of the full theory given in Section 2.1.2, where we found that the 
/3 = 0 theory has four free parameters, in order to see that the R = 0 condition imposes 
two parametric constraints on generic static spherically symmetric solutions. 

The R = 0 solution for the (0,0) family is obtained if and only if b 2 = a 2 - 


A{r) = l + a 2 (r 2 + r 


4 12 aa 2 + 7 6 320a 2 a2 + 100 aa 27 + 7 


B(r) 

b 0 


= 1 + a 2 r 2 + r 


20a 

~ 4 24aa 2 + 7 
40a 


+ r 


+ r 


1120 a 2 
6 960a 2 a 2 + 144 aa 27 + 7 2 


3360a 2 


+ 0(r 8 ) 
+ 0(r 8 ) 


(6.9) 

( 6 . 10 ) 


and R^ y = 0 is obtained if and only if 02 = 0. Since A(r —> 0) —> 1 for all of (0, 0) this 
solution must have an even number of horizons. 

The R = 0 solution for the (1,-1) family is obtained if and only if 04 = ((0162 — of: 


ajb 2 J+0(r 6 ) (6.11) 

= - + ai + b 2 r 2 + ^a\b 2 r 3 - \r A a\b 2 + 0(r 5 ) . ( 6 . 12 ) 

6 _i r 65 

The free parameter a\ corresponds to the Schwarzschild mass and the free parameter b 2 con¬ 
trols the deviation from the Schwarzschild solution. Specifically, one has pure Schwarzschild, 
i.e. R^ u = 0 , if and only if b 2 = 0 , and inspection of the solution then shows that 
01 = — ^(Mschwarzschiid) -1 ; so we expect the solution to have a horizon for ai < 0 and 
to have no horizon for <27 > 0 . 

The (2, 2) family (which does not appear in general relativity) has a four-parameter 
family of R = 0 solutions but cannot have R^y = 0 for all r, because, e.g., in the (2,2) 
family one has R rr = 3r -2 + ©(r^ 1 ). 

A summary of boundary structure and parameter counts for the three near-origin solu¬ 
tion families is given in Table [2j 


A(r) = a±r — a\r 2 + a 3 r 3 + r 4 ( ^-a\b 2 — af ) + r' 


, 3^3 


a i 
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Table 2: Trace no-hair boundary structure and solution parameter counts 
(including the trivial time-rescaling parameter) 


(s,t) solution family 

VABr 2 Rd r R 

number of free parameters 

generic theory 

/3 = 0 theory 

( 0 , 0 ) 

0 (r 3 ) 

3 

2 

(1,-D 

0 (r 3 ) 

4 

3 

( 2 , 2 ) 

Oir- 1 ) 

6 

4 


7 Expansion around a nonzero radius ro 


We can refine our understanding of the various forms of solutions by expanding now around 
an arbitrary radius tq. This will be easier to do in a slightly different set of variables: 


ds 2 = — B(r) dt 2 + + r 2 d£l\ 

f(r) 


(7.1) 


related to the usual Schwarzschild variables (1.3) by A(r) = l//(r). 


We can use a Frobenius ansatz for the expansion about ro similar to our expansion (3.1) 
about r = 0 : 


/ =fw(r - r 0 ) w + f w +i(r - r 0 ) 


B 


=(r - roY + b t+ i{r - r 0 ) t+1 + . 


r +l +... 

(7.2a) 

+ ... 

(7.2b) 


for some exponents w and t, not confusing these undetermined (w, t) exponents with the 


undetermined (s,t) exponents used earlier in the expansion (3.1) around r = 0. We shall 
find various Frobenius solution families and also some non-Frobenius families of solutions. 
However, detailed discussion of all of these is beyond the scope of this paper. In Section 


7.2 we shall present a summary and discussion of all the solutions that we have found. For 


the purposes of the main thread of our discussion the situation where the space-time has a 
horizon will be of most interest to us and we turn to it now. 


7.1 Solution family around a horizon 

We now focus on the properties of spherically-symmetric solutions with horizons. We saw in 
Section [ 6 ] that in a static asymptotically flat spacetime the presence of a horizon implies that 
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the Ricci scalar must vanish and consequently the system becomes equivalent to Einstein- 
Weyl gravity, with Lagrangian 


e~ l C = ^R- aC pvprT C, 


/iispa • 


(7.3) 


Using the metric (7.1) the equations of motion then imply 


0 =r 2 (B (2 fB" + B'f ) - fB ,2 ) 

+ 4 B ( rfB 1 + B (rf + / - l)) (7.4a) 

0 = a ^ r 3 fB' l-BB’f + fB 12 - 2 B 2 f"] 

- 2 B 2 f ( r 2 B’f - 2 B (r 2 /" -rf + 2)) - Bf 2 (3 r 2 B' 2 + SB 2 ) 

- rB 3 f (3 rf - 4) j + 2 7 r 2 R 2 (rfB' + B(f - 1)) . (7.4b) 


The Schwarzschild metric itself is of course still a solution to Equations (7.4), with 


B(r) = f(r ) = 1 - — . 

r 


(7.5) 


In the higher-derivative theory we do not have a general solution in closed form but a 
Frobenius analysis performed around the horizon can reveal its relation to the Schwarzschild 
solution. We consider solutions in the neighbourhood of a horizon, assumed to be at r = r o- 
By definition, the metric function B(r ) vanishes at a horizon so we look for solutions of the 
fompl ( |7.2p with t > 0. Using the expansions of B and / as given in ( |7.4b| ) shows that 
w < |. In ( |7.4a ), for w < | it is the terms in the first line that contribute the leading-order 
term in (r — ro), and the reader can easily check that these vanish only if t = 2 — w. Using 


t = 2 — w in (7.4b) then shows that the first line contributes the leading order term in 


(r — r o) and this then vanishes only if w = t = 1. Thus both B and / must be linear in 
(r — ro) at the horizon, as in the Schwarzschild case. 

The solution to the equations of motion then is found to have three free parameters: 
6 i,/i,r 0 . The other coefficients b n , f m can be solved for in terms of these: 


(7.6a) 


(I) - 8^ + ^ + 4) (r-r 0 f + O((r-r 0 r), 

{~^T + o 72 + 7“2 “ ~) ( r “ r o) 2 + 0((r- r 0 ) 3 ) , (7.6b) 

V 8a fi 8a ffr 0 f\rf r 0 J 


f =Mr-r o)+(g-^ 

f =(r - r0) + (“to/, 


2 Note that an expansion of the form for an asymptotically flat solution with a horizon, i.e. a 

(§’ 2 )y r - r expansion in the notation of Table |3j is not possible, as one must have vanishing Ricci scalar 
for asymptotically flat solutions and such an expansion does not then exist. 
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and so on. The parameter b\ is trivial, as we have seen, in the sense that it can be absorbed 
into a rescaling of the time coordinate. Note that the Schwarzschild solution corresponds 
to the case where f\ = 1 /ro, with the same time coordinate as for Schwarzschild if one also 
chooses b\ = 1/ro- 

This count of the free parameters tells us which of the three families near the origin 
corresponds to asymptotically flat solutions with horizons. Consider the (s,t) = (1,-1) 
family of Table [2} As we have seen in Section [6j for asymptotically-flat solutions this family 
must have R = 0 and its equations of motion accordingly must become equivalent to the 
(5 = 0 theory as considered in this section. In the (5 = 0 theory, one accordingly has three free 
parameters as shown in Table [2] The theory thus contains the two-parameter Schwarzschild 
solution, together with a one-parameter family of deviations from Schwarzschild. Nearby 
the Schwarzschild solution within this one-parameter family of deviated solutions, we would 
certainly expect the horizon still to be present. So we expect the (1, —1) family to be a 
three-parameter family, in which a three-dimensional volume of that parameter space has 
a horizon. Therefore in the (5 = 0 theory we identify the ( s,t ) = (1, —1) indicial solution 
family obtained from expansion near the origin with the solution family containing horizons. 

One needs to be careful with the logic here. Although we say that asymptotic flatness 
implies R = 0 for the (s,t) = (1,-1) solution family, and therefore its equivalence to 
solutions of the Einstein-Weyl theory, the converse is not necessarily true. It is likely that 
not all (1,-1) solutions of Einstein-Weyl gravity are asymptotically flat, and in fact the 
loss of asymptotic flatness is a natural guess for the consequence of turning on the ‘non- 
Schwarzschild’ parameter measuring the deviation from Schwarzschild. In Section [8j we 
will further explore such deviations from general relativity in Einstein-Weyl gravity and 
the implications for their asymptotic behaviour through a linearised expansion in the non- 
Schwarzschild parameter. 


7.2 Summary of expansion behaviours around a non-zero radius ro 


The Frobenius ansatz (7.2) has three solution families. The first is simply the (0,0) family, 
corresponding to no special radius. Its main value is reinforcing our conclusions about the 
generic number of free parameters of the theory. The second is the (1,1) family already 
discussed, describing a horizon. The third is the (1,0) family, which we shall describe as a 
wormhole in Section [9l 

As mentioned in our Section [3] discussion of expansions around the origin, a natural 
concern is that there might be other solution behaviours not captured by the integral-step 
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Frobenius expansion. We cannot check all alternative expansions exhaustively, but we have 
tried a variety of non-Frobenius expansions and have found two other solution families. Both 
of these non-Frobenius families involve half-integer as well as integer powers of (r — ro). As 
before, we denote Frobenius families by bracketed pairs of indices (w, t), and we shall denote 
the new series similarly but with a subscript: (w, where they go as 

/ =/o(r - r 0 ) w + /i(r - r 0 ) w+ 5 + f 2 (r - r 0 ) w+1 + ... (7.7a) 

=(r - r o y + bi(r - r 0 ) t+ 5 + b 2 {r - r 0 ) t+1 + ... (7.7b) 

t>o 

We shall not give a detailed analysis of all five solution families in this paper, but we 
summarise our findings in in Table [3| 

Table 3: General expansion behaviours around a nonzero radius ro 


(w,t) solution family 

\J~AB r 2 Rd r R 

number of free parameters 

generic theory 

j3 = 0 theory 

(0,0) 

0 (1) 

6 

4 

(1,1) 

0(r - r 0 ) 

4 

3 

(1,0) 

0{y/r - r 0 ) 

3 

2 

(1,0) 

0 (1) 

6 

4 

(3 1\ 

V 2 ’ 2 / yjr—ro 

o 

1 

*1 

o 

3 

N/A 


Note that the (|, |) solution family does not occur in the /3 = 0 theory. The (1,0) 
solution family is a subset of the (l,0)^rr^ family, obtained by setting the coefficient of 
odd powers of \Jr — ro to zero while holding the coefficients of even powers at finite values. 

8 No-hair theorem for a linearised deviation from Schwarzschild 

We saw in Section [6] that by considering the trace of the field equations for gravity with 
general quadratic curvature terms added, one can derive a no-hair theorem that shows that 
the Ricci scalar must vanish in any asymptotically-flat black hole solution. Unfortunately, 
similar arguments applied to the full set of field equations fail to establish a more powerful 
result that one might have hoped to demonstrate, namely the vanishing of the full Ricci 
tensor for all asymptotically-flat spherically-symmetric solutions with horizons. Had one 
been able to obtain such a result, this would have shown the Schwarzschild solution to be 
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the unique static spherically-symmetric asymptotically-flat black-hole solution in theories 
of gravity with curvature-squared corrections. 

Indeed, as we found in Ref. m , there are non-Schwarzschild black-hole solutions to be 
found numerically, so the failure of a full Lichnerowicz-type no-hair theorem including the 


traceless components of the field equations ( 2 . 1 ) is now seen to have been quite indicative. 
Nonetheless, one may still obtain useful information about the set of solutions with horizons 
from a no-hair theorem analysis carried out to linearised order in the non-Schwarzschild 
parameter discussed in the last section. As we shall see, the upshot from this analysis is 
that, provided the curvature-squared terms have sufficiently small coefficients in comparison 
to the scale size of the black hole, then there can be no well-behaved static and spherically- 
symmetric black holes that are perturbatively close to the Schwarzschild solution. So in 
this restricted sense, one can show that the Schwarzschild solution is generally an isolated 
solution, discretely separated from other asymptotically-flat solutions with horizons. 


We may consider solutions of the equations (7.4) that are infinitesimally close to Schwarzschild 
by writing 


B(r) = 1 — 




1 


= f(r) = 1 - 


r o 


( 8 . 1 ) 


r(l + eZ B {r)) A(r ) J v ' " r(l + eZ A [r)) ’ 

and keeping only terms of order e. From the two coupled equations of motion in Z A and 
Zb, a single 3rd-order ordinary differential equation purely for Z A {r) can be obtained by 
eliminating Z B (r): 


Z B (r) - (r - r 0 )Z B (r) = Z A (r) + 

+ 


a (—8r 2 + 16rr 0 - 9 r%) ( r - r 0 ) ^ 

2 yr 4 — 2yr 3 ro — 4arro + 5 otq 

2ar(2r-3r 0 )(r-r 0 ) 2 „, ^ 

- ■ 


( 8 . 2 ) 


2 yr 4 — 2yr 3 ro — 4arro + 5 ar^ 

In fact, the resulting equation in Z A involves only Z' A , Z" A and Z" A terms, and consequently 
we have a 2nd-order ordinary differential equation for Z' A . It is useful to introduce a new 
variable Y(r), defined by 


z A{r)= / Y(r)u(r)dr, 


(8.3) 


'r 0 


where u>(r) = 1 for now but this will be revised later. The lower limit in (8.3) is chosen to 
be ro in order to ensure that Z A {r) vanishes on the horizon 13 Using the abbreviation 


C = a ('T r o) \ (8-4) 

1,i The third linearly independent solution of the 3rd-order ordinary differential equation for Za{ r) itself is 
Za (r) = const., which, as can be seen from ( |8.1| ), just describes a perturbation that shifts the mass of the 
Schwarzschild solution by an infinitesimal constant. 
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the second-order ordinary differential equation for Y (r) is then 


h 0 Y + hi Y' + h 2 Y" = 0 , 


(8.5) 


where one has 


ho = 2r 7 — 2ror 6 — 8rgr 5 C + 16rQr 4 £ — 5 r$r 3 ( — 32rgr 2 £ 2 + 44rgr<C 2 — 20 r 7 0 C, 2 , 

hi = 4r 2 (2r - 3r 0 ) r%( (r 3 - r 0 r 2 + rgC) , 

h 2 = — 2r 2 (r - r 0 ) ?’gC ( 2 (?’ - r o) f’ 3 + ?’o ( 5r o - 4?’) C) • 


( 8 . 6 ) 


One can easily see from (8.5) that at large r the two solutions to the field equations go 


like 

Y ~ a\ (1 + m 2 r) e _T7l2r + a 2 (1 — m 2 r ) e m2r (8.7) 

where m 2 = \J r y/2a as before. For r close to the horizon at r = ro, one can once again 
use the Frobenius method to find the r —> ro asymptotic behaviour of the two independent 
solutions Y\ and Y 2 . We find that they take the asymptotic forms 


Y\ = 1 + ci (r — ro) + c 2 (r - ro ) 2 + C 3 (r - ro)" 

b. 1 


Ci = Ci(a, r 0 ) for i > 1 


Y 2 = Yi log(r - r 0 ) H-— + h (r - r 0 ) + b 2 (r - r 0 ) 2 + b 3 (r - r 0 ) 3 H-, 

r - r 0 
a 

with 6 _i =-, = bi(a, ro) for i > 1. ( 8 . 8 ) 

ro 

Thus, in order for the metric perturbation ZA{r) to be non-singular at large r we must 
have cl 2 = 0, while for non-singular behaviour near r = ro the overall coefficient of the 
Y 2 solution must be zero. We shall now attempt to show that no such solution Y that 
interpolates between these limiting forms can exist. 


To do this, we take the Y equation (8.5), multiply it by u{r)Y{r) for some chosen u{r) 


and then integrate the result from r = ro (the horizon) out to infinity. Firstly, we note that 
0 = (ho Y+hi Y'+h 2 Y")uY = uhoY 2 -uh 2 Y' 2 + (uh l -u'h2-uh' 2 ) YY'+(uh 2 YY ')'. (8.9) 


We now choose u(r) so that the coefficient of YY' vanishes, by solving 

uh\ — uh 2 — uh' 2 = 0 . 

This gives, up to a constant factor that we may take to be 1, 

(r - ro) 


u(r) = 


(2 (r — ro) r 4 + r[j (5ro — 4r) r£) 


2 • 


( 8 . 10 ) 


( 8 . 11 ) 


32 







Integrating (8.9) from ro out to infinity, we thus obtain 


dr(uh 0 Y 2 - uh 2 Y'~) = 0. 


( 8 . 12 ) 


'ro 


Note that we can drop the total derivative term, since we have established that an acceptable 
solution for Y must satisfy Y ~ const, near r = ro and Y ~ e~ mr near infinity. If we can 
show that uho and —uh 2 are both non-negative in the interval ro < r < 00 , then we will 
have shown that no acceptable solution Y can exist. 


The function u{r ) obtained in (8.11) is manifestly non-negative in the range tq < r < 00 . 


It is then evident from ( |8.6[ ) that showing the non-negativity of uho and —uh ,2 is equivalent 
to showing that ho and H 2 , given by 


ho = 2 r‘ — 2ror 6 — 8r^r 5 ( + 16rdr 4 £ — 5 r^r 3 ( — 32rdr 2 £ 2 + 44r®rC 2 — 20rg(( 2 , 

H -2 = h2 —r 27 = 2 (?’ — r 0 ) r 3 + rg (5r 0 — 4r) ^ (8.13) 

-2 r 2 ( r - r 0 ) rgC 

are non-negative in the interval rg < r < 00 , for some range of ( > 0. 

It is easy to see that H 2 is non-negative in the whole interval if and only if 

0 < C < y • (8-14) 

The non-negativity of ho provides a stronger condition on £. Setting r = rg, we see that 
h 0 (r 0 ) = rg (3 — 8 C) Cj an d so we must certainly have 0 < C < §■ An easy way to investigate 
the bound on ( under which ho(r) is non-negative in the entire range rg < r < 00 is to write 

r = r 0 (l + x), £ = ( 8 . 15 ) 

1+2/ 

since then the ranges ro < r < 00 and 0 < C < Cmax are mapped into the positive quadrant 
0<x<oo, 0<y<oo. We then bring the expression for ho over a common denominator 
(which is the manifestly positive quantity (1 + y) 2 ), and examine the numerator, which is 
a multinomial in x and y i.e. it is of the form 

N M 

Y J Y, c ^ xn y m ( 8 - 16 ) 

n m 

where the C n ^ m are functions of ro and Cmax- The condition C n ^ m > 0 for all n, m is clearly 
sufficient (but may not be necessary) for non-negativity of ho in ro < r < 00 . This condition 
easily yields the bound 

0 < C < ^ • (8-17) 

Conversely, as we have seen, if £ exceeds this bound then ho will be negative at r = rg. 
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The upshot from this discussion is that, provided £ is bounded from above by (8.17), 


then there cannot exist a regular infinitesimal perturbation of the metric away from the 
Schwarzschild solution. For £ exceeding this bound, we can gain no information from this 
discussion. 

One can actually improve the upper bound on ( by making a different choice of ui(r) 


in (8.3), and by repeating the previous steps. For example, if we take u(r) = (cro + r) _1 , 
where c > — 1 is a constant to be chosen, the optimal bound £ max is the largest positive 
root of the sextic 

9600 CLx + 8624 cLx - 9360 dx - 4461 (Lx + 1216 Cmax + 1116 Cmax - 48 , (8.18) 

which is approximately given by Cmax ~ 0.617292. This is achieved by choosing the constant 
c to be given by 

C = ~ 36 Cmax + 19 Cmax + 2 ^ 0 , 164789 . (8.19) 

20 Cmax — 9 Cmax — 2 

A slight improvement on this can be achieved by taking instead u(r) = ( cr jj + r 3 )-i /3 
We then find 0 < C < Cmax with Cmax the largest positive root of 


which is Cn 


28160 Cmax + 12176 Cmax - 43374 Cmax + 19179 Cmax ~ 2322 = 0 , 

~ 0.6262615, attained when the constant c is chosen to be given by 
3 - 4 Cma 


c = 


1 Cmax 3 


0.2462346. 


( 8 . 20 ) 


( 8 . 21 ) 


The condition of C n m > 0 in (8.16) is a sufficient condition, and the necessary bound 


on C may be better. Trying different functions cj(r) could presumably improve the bound 


further. The best one could hope to achieve by this method is the bound (8.14) (valid for 


any function w(r)) arising from the need for H 2 to be non-negative also. In any case, we 
have established that, provided C is sufficiently small, there is a linearised no-hair theorem 
that rules out regular black holes that are close to the Schwarzschild solution. 

The treatment of the traceless components of the field equation by such a linearised 
perturbative approach establishes, for C appropriately bounded, that the Schwarzschild 
solution is isolated : there can be no other nearby asymptotically-flat solutions with horizons. 
The notion of ‘nearby’ solutions is made clear by our general analysis of the parametric 
dependence of solutions with horizons. From the parameter count summarised in Section 
??, we presented an argument in Section[7]for identifying the 3-parameter family of solutions 
with horizons with the 3-parameter (1, —1) family of solutions around the origin. The classic 
2-parameter Schwarzschild solution is clearly contained in this family. It is in terms of the 
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single remaining parameter of the ( 1 , — 1 ) family that one can characterise the ‘distance’ of 
other nearby solutions with horizons from the Schwarzschild solution. 

The need to make such a cautious statement about the isolation of the Schwarzschild so¬ 
lution, however, clearly raises questions as to whether there might exist other asymptotically- 
flat solutions with horizons that are not ‘near’ to Schwarzschild. From Ref. m , we know 
this in fact to be the case. In the next section, we give a brief overview of what can be said 
about such solutions from numerical studies. 


9 Numerical analysis 

The detailed nonlinear field equations (given in AppendixjA]) for our Schwarzschild-coordinate 
spherically-symmetric system in the general a, /?, 7 theory ( 1 . 1 ) are clearly not very 
amenable to a closed-form solution. Having studied the asymptotic behaviour of solu¬ 
tions at the origin, at spatial infinity and at a horizon in Sections m and[7j we now need 
to consider what happens in between these various limiting regions. This is only approach¬ 
able by numerical study. We do not purport to give an exhaustive treatment of numerical 


solutions to the theory ( 1 . 1 ) here, but some review of what is already known and what can 
be obtained by Mathematica experimentation is in order m- 


9.1 (2,2) solutions 

Firstly, let us consider solutions that could be obtained from coupling to a positive-energy 
shell source as discussed in Section [5j Only the (2, 2) indicial family has the full count of 
six parameters that are required to satisfy the six continuity and jump conditions across 
a shell source. After such matching, two parameters must implicitly be used to guarantee 
the absence of rising exponential behaviour at spatial infinity, corresponding to the rising 
spin-two and spin-zero Yukawa terms of the linearised theory. It is not known, however, 
which combinations of free parameters near the origin, given in Table [lj need to be tuned 
so as to eliminate the rising behaviour at spatial infinity. In order to match a (2, 2) family 
solution on to asymptotically-flat behaviour at infinity, one procedure is to start with a 
series-expanded solution near the origin and integrate outwards numerically, and also to 
start from an asymptotically-flat solution at spatial infinity and integrate inwards, then 
to adjust parameters so as to make the two numerical solutions match at an intermediate 
radius. Such a procedure was carried out in Ref. HD for the theory with m 2 = mo, which in 
the notation of this paper means a = 3/?. Owing to the trace no-hair theorem as presented 
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in Section [gJ any asymptotically-flat solution that has any amount of falling spin-zero 
Yukawa behaviour near spatial infinity cannot have a horizon; as one can see from Eq. 
(4.8), such solutions necessarily have R ^ 0 in the r —> oo asymptotic region. Indeed, 


the asymptotically flat (2,2) family solution found in Ref. |llj displays a dominant 1/r 
Schwarzschild-type behaviour as r — > oo, but it also displays a falling Yukawa correction 
and deviates strongly from Schwarzschild at smaller r. It does not encounter a horizon at 
intermediate r values, but limits to (2, 2) family behaviour near the origin. The calculation of 
Ref. DU was made for a normal positive-sign mass M (given [20] by 87r times the coefficient 
of 1 /r in gtt as r —> oo for a theory with 7 = 1 ). 

A similar calculation can be made in the 7 R — aC 2 theory with f3 = 0, in which a 
vanishing Ricci scalar, R = 0, is naturally guaranteed. Accordingly, this theory benefits 
from a further reduction of the set of third-order nonlinear field equations as given in 
Appendix [A] down to a pair of second-order equations. These are equivalent to the system 
( |7.4[ ) for the f = 1/A and B variables. This f3 = 0 theory also has asymptotically flat, 
limiting to (2,2) indicial, solutions without a horizon. A representative numerical solution 
with positive mass M is shown in Figure [l] 



-InA 



Figure 1: Horizonless asymptotically-flat solution at spatial infinity limiting to a (2,2) 
indicial family solution at the origin. The right-hand plot for — lnA(r) and 
lnR(r) shows the relation to the Schwarzschild solution. 


The numerically obtained asymptotically flat (2,2) solution shown in Figure [I] was ob¬ 
tained for a positive ADM mass M but with a negative C 2 - falling asymptotic Yukawa 
coefficient. Comparing to linearised theory solutions obtained with various sorts of distri¬ 
butional matter coupling as discussed in Section [4j one sees that this relative sign between 


M and C 2 - differs from that obtained for a point source in (4.7). Recall, however, that 


for the higher-derivative theory (1.1) there is no Birkhoff theorem, and one does not have 


36 















any expectation of a universal sign relation between M and 6*2-. Comparing instead to 
the linearised theory with the ‘balloon’ source ( 4 . 18 | ) , one finds that 62- can have either 
sign compared to M, so the sign found for C2- in the Figure [l] solution poses no particular 
puzzle. 

While in Ref. m the method was to combine a shooting-out calculation from the origin 
with a shooting-in calculation from large radius, in order to find solutions like that of Figure 
[l]we employed a simpler method of just integrating inwards from a large radius using initial 
conditions taken from an asymptotically flat solution of the linearised theory. In general, 
this leads to divergent behaviour for A{r ) or B(r ) at small radii, but one should take care to 
notice that there are two kinds of divergent behaviour that can occur. Holding the linearised- 


theory’s coefficient 62,0 ~ —M fixed in ( 4 . 2 ) while varying the falling Yukawa coefficient 
C2-, one finds ranges of C2- values for which A(r) —> 00 while B(r) —> 0 as r —> 0 , but 
then one finds a different range of C2- values for which A(r) —> 0 while B(r) —> 00. In 
between these ranges one finds (subject to numerical accuracy) a value of C2- for which 
both A(r) —> 0 and B(r) —>• 0 : this is a ( 2 , 2 ) solution as shown in Figure [lj 

This procedure for finding asymptotically flat ( 2 , 2 ) solutions such as that of Figure [I] 
reveals another feature of the overall solution space that calls for further study. If the ( 2 , 2 ) 
indicial family solution lies on a separatrix in parameter space between two other kinds of 
more generic solution, what are these other kinds of solution? One of these other kinds of 
solution may be viewed as a ‘wormhole’, to which we turn our attention next. 


9.2 Wormholes 

Another type of solution that can be found numerically may be described as a ‘wormhole ■0 
Such solutions are characterised by the existence of a zero for f(r) = 1 /A(r) but with 
B(r) = —gtt 7^ 0 . We have seen in Section [ 7 ] that solutions with B(r) vanishing at some 
radius ro must also have f(r) vanishing at r0 as well. However, the converse is not necessarily 
true. Integrating inwards from an asymptotically flat solution at spatial infinity, one finds 


such solutions starting from a linearised solution ( 4 . 2 ) with chosen values of 6*2,0 and 62 _. 
In this way, one finds solutions with either sign of M = —8716*2,0 and either sign of the large 
r spin-two Yukawa coefficient 6*2 - • 

Another way to investigate such solutions with /(ro) = 0 but B(tq) / 0 is again to use 
Frobenius asymptotic analysis to find the possible behaviour as r —> ro and then to integrate 
outwards, looking for asymptotically flat solutions. Asymptotic analysis indeed shows, as 


1 Wormholes have also recently been considered in the pure R 2 theory in Ref. EH. 
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one can see from the existence of the (w, t ) = (1, 0) and (w, t ) = (1,0)^rr^ solution families 
shown in Table [3j that there can be solutions for which f(r) vanishes at some radius ro 
but where B(r) remains at some nonzero value. We have studied this in particular in the 
7 R — aC 2 theory with (5 = 0. 

As we have seen above in Section [7.2| and summarised in Table [3j asymptotic analysis as 
r — > ro for f(r o) = 0 but with B(tq) 0 turns up the following situation. The leading term 
in fir) is always linear in (r — ro) and the leading term in B(r ) is always, by assumption, 
a constant. Since we are also by assumption considering the first zero of /(r) as r comes in 
from infinity, without B(r ) having yet crossed zero (which would have constituted a horizon 
as we saw in Section [ 7 ]), the B(ro) = bo constant must be positive. It is at this point that 
the half-integral-step expansions of type ( w,t .) = (l,0) v ^rr^ from Table [ 5 ] become relevant. 
If half-integral steps are allowed, then one finds an expansion with four free parameters, 
which is the generic number of free parameters for spherically symmetric solutions in the 
7 R — aC 2 theory, which has two second-order field equations (equivalent to those given in 


15 


(7.4)) for f(r) and B(r) 

Assuming instead integral steps in powers of (r — ro) after an initial zero at ro leads 
to a more constrained solution system with only two free parameters: the trivial rescaling 
parameter affecting B = —gtt (which may be realised as bo) and ro itself. Numerically 
integrating outwards in r, one then generally finds rapidly divergent behaviour as r —> oo, 
but this behaviour can be of two different types, similarly to the two types of divergent 
behaviour surrounding the asymptotically flat (2,2) solutions discussed in Section 9.1 one 
type has /(r) —> oo and B(r ) — > 0, and a different one has /(r) — > 0 and B(r ) —> oo. In 
between these behaviours, by tuning ro one can find a solution that becomes asymptotically 
flat for a specific value ro = r*. Such a solution is shown in Figure [2} Comparing the r —>• oo 
asymptotic behaviour obtained by numerical calculation to the asymptotically-flat case of 


the linearised theory solution (4.2) with (5 = 0, one finds such an integral-step solution 
corresponding to M 6 * 2,0 < 0 and C* 2 _ < 0 . 


15 A confirmation of this analysis may be observed in numerical solutions shooting inwards from asymp¬ 
totically flat spacetime. For /(ro) = 0 but B(ro) > 0, there is an apparent divergence in the gra¬ 
dient of B as one approaches ro, where / has a zero, agreeing with an expansion structure /(r) = 
/i(r - r 0 ) + / 3 / 2 (r - r 0 ) 3/2 + ... and B(r) = b 0 + b 1/2 {r - r 0 ) 1/2 + .... 
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Figure 2: Z2 symmetric wormhole solution with f(r) —> 0 as r —> ro (dashed line) and 
B(r ) —> const, (solid line). The left plot shows the small-scale structure near 
ro — 0.57 and the right plot shows larger-scale structure. 


To see why such a solution may be described as a wormhole, consider /(?’o) = 0 and 
B{tq) = bo > 0 and expand in (r — ro): 

f(r) = (r - r 0 )/'(r 0 ) 4-, B(r) = b 0 + B'(r 0 )(r — r 0 ) 4-. 

As one can see from the left-hand calculated plot in Fig. [2j both B'(ro ) and /'(ro) are 
positive. Now make a coordinate change in the radial coordinate according to 

r-r 0 = \p 2 , (9.1) 


after which one has the asymptotic form of the metric 

ds 2 = —(bo + \B'(ro)p 2 )dt 2 + + (rg + \rop 2 )dQ 2 + ... . (9.2) 

Since this solution is an even function of p, it is naturally Z2 symmetric in p and can be 
extended to the full range —00 < p < +00. Geodesics in the p > 0 patch match smoothly 
onto geodesics for p < 0 and so continue on naturally through to negative p without hitting 
a singularity. 

The interpretation as wormholes of the more general non-Z2-symmetric solutions with 
/(r 0 ) = 0, B(ro) = const, arising from asymptotic expansions with half-integral steps in 
(r — r 0 ), i.e. in the (w, t ) = (1,0)^rr^ family, is less clear. As we have seen, such solutions 
have four parameters in the expansion about ro, which is the generic number for spherically 
symmetric solutions in the -yR — aC 2 theory. This could allow tuning of one parameter 
combination so as to ensure asymptotic flatness at spatial infinity, even for an arbitrary 
value of ro- The expansion in half-integral powers of (r — ro), however, leads to odd powers 


of p after making the coordinate change (9.1). This destroys the Z2 symmetry of the ?’o- 
tuned solution and invites the question whether one will then have B —> 0, and consequently 
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a horizon, at some value rhor of the radius. Accordingly, the interpretation of such general 
Z 2 -asymmetric solutions as wormholes is not so clear as for the Z 2 symmetric solutions. 


9.3 Schwarzschild and non-Schwarzschild black holes 


Turning now to asymptotically-flat solutions including a horizon, we know from the trace 
equation no-hair theorem of Section [ 6 ] that all such solutions must have vanishing Ricci 
scalar, R = 0. Accordingly, the analysis of such solutions can be restricted to the jR — aC 2 
theory, since the field equations of the general 7 , a, (3 theory 0 reduce to those of the 
/3 = 0 theory when R = 0. Furthermore, the results of Sections [ 7 ] and [ 8 ] show that the 
Schwarzschild solution is in general isolated in the sense that the linearised no-hair theorem 
of Section [ 8 ] does not permit, for 0 < ( = a( 7 rg ) -1 < £ max , where ro is the horizon radius, 
any solution infinitesimally deviating from Schwarzschild in the single non-Schwarzschild 
parameter allowed by parametric analysis near the horizon, as explained in Section [7j As 
presented in Ref. [15j . however, the qualified nature of the linearised no-hair theorem led to 
a suspicion that there might in fact be other asymptotically-flat solutions with horizons that 
are in general distinctly separated from Schwarzschild in the value of their non-Schwarzschild 
parameter. In the notation of Ref. [Tf)!, the non-Schwarzschild parameter may be taken to 
be 5, defined by 


f( r ) 

h 


vtt = fi (r - r 0 ) + 0(r - r 0 ) 2 
A(r) 

1 + 5 

ro 


(9.3) 


The limiting value Cmax of the linearised no-hair theorem suggests the existence of a branch 
point in the value of the horizon radius ro = r™ in at which infinitesimal values of <5 can 
give rise to non-Schwarzschild asymptotically-flat solutions with horizons. As found in Ref. 
[15], for ro > r™ in such non-Schwarzschild solutions do exist, but they lie outside the linear 
validity range of the 5 parameter expansion. So, except at the r“ m branch point, the 
Schwarzschild solution must be isolated. 

Numerical calculations reveal the properties of the various black-hole phases. The phase 
structure m in terms of black-hole mass GM (as above, given [20j by | the coefficient of 
1/r in gtt as r —> 00 ) is shown in Figure |3j 

Joining the non-Schwarzschild r > tq solution outside the horizon to the r < interior 
solution, one obtains the result shown in Figure [4] The deviation from a Schwarzschild 
black hole is shown by the fact that f(r) = 1 /A{r) ^ 5(r); however, the approach as 
r —> 0 shows this solution can still fit naturally into the (1,-1) indicial family. Another 
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Figure 3: Phase structure of the Schwarzschild (dashed line) and non-Schwarzschild (solid 
line) black holes in a theory with a = sketched for a theory with 
G = = 1. The Schwarzschild mass is given by GM = fro- 

(B(r),f(r)) 



Figure 4: Non-Schwarzschild black hole for GM ~ 167r 0.276 with a horizon at r = 1. The 
dashed line denotes B(r) and the solid line denotes f(r) = 1 /A{r). 


piece of information obtained from the calculation producing Figure [4] is the sign of the 
C 2 - coefficient of the falling e ™ 2r spin-two Yukawa term. Both for positive and negative 
GM = — \C 2 fi, the sign of C2- appears to be negative for such solutions. For M > 0, this 
sign is opposite to that which would have been expected from the linearised 7 R—aC 2 theory 


coupled to a positive-energy shell delta-function source, as shown in (4.7a) or (4.12a). 

The numerical results presented in this section are clearly only an initial foray into the 


perhaps rich phase structure of the solution space of theories derived from the action (1.1), 
and this subject clearly requires more careful numerical analysis. 


10 Conclusion 


In this paper, we have carried out an analysis of the static spherically-symmetric solutions 


of the field equations derived from the action (1.1). This extends older results [3j by a full 
asymptotic analysis of the indicial (s,t) = (0,0), (1,-1) and (2,2) solution families near 
the origin, together with a careful count of the parameters occurring in each family. The 
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difficult question is then what happens in the intermediate 0 <r<oo region before one 
reaches spatial infinity, near which the assumption of asymptotic flatness once again allows 
for a closed-form study of solutions via the linearised version of the field equations derived 


from (1.1). 

In Section [6j we reviewed the no-hair theorem for the trace of the higher-derivative 
theory’s field equations, agreeing with this part of the results of [13]. This implies that 
asymptotically flat solutions containing a horizon must have R = 0 throughout the whole 
extra-horizon spacetime. Consequently, one knows that asymptotically-flat solutions dis¬ 
playing a spin-zero Yukawa term ( i.e. with a non-zero coefficient Co-) cannot have a hori¬ 
zon, since Co- ^ 0 ==> R ^ 0. Study of solutions with vanishing Ricci scalar R, and 
hence without such a spin-zero Yukawa term, reduces to the somewhat simpler case of the 
7 R — aC 2 theory with (3 = 0. In this restricted context, the field equations can be reduced 
to a system of two coupled equations with at most second derivatives in ^. Although there 
is no general no-hair theorem for the trace-free components of the field equations, one still 
does have a linearised no-hair theorem for the traceless field equations as derived in Section 
[8] This implies that the Schwarzschild solution is generally isolated within the family of 
(1, —1) solutions, which is controlled by a single non-Schwarzschild parameter, provided the 
horizon radius ro is larger than a certain bound y/ot/( yCmax)- As one moves away infinitesi¬ 
mally from the Schwarzschild solution within the (s,t) = (1,-1) family, the only thing that 
can generally happen for solutions with a horizon is that asymptotic flatness is lost. 

What happens outside the domain of validity of the linearised no-hair analysis is an¬ 
other matter, however, and in Ref. m asymptotically-flat non-Schwarzschild black-hole 
solutions were indeed found. Except near the ro > \/«/( 7Cmax) horizon radius bound for 
the linearised no-hair theorem, such non-Schwarzschild black holes can only exist owing to 
nonlinear dependence on the (1,-1) family non-Schwarzschild parameter. Numerical evi¬ 
dence points to the ro ~ V^TCrCmaxJ horizon radius coinciding with a branch point in the 
black-hole solution space, at which the non-Schwarzschild black holes first occur (c./. Figure 
[3]). This is clearly consistent with a breakdown of the linearised no-hair theorem for that 
radius, because just above such a branch point the non-Schwarzschild black holes should 
indeed be obtainable by linearised perturbation away from the Schwarzschild solution. 


The general structure of the solution space for the higher-derivative (1.1) theory is still 
not completely clear, but we have found evidence of a rather rich phase structure for the 
static spherically symmetric and asymptotically flat solutions of this theory. We have seen 
from analysis of shell continuity and jump conditions that only the (s,t) = (2,2) solution 
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family (without horizons) of Section |9.1| can couple to a standard distributional matter 
shell, but this need not imply that the (2, 2) family is the only equilibrium endpoint of 
gravitational collapse. In addition to the (2,2) family, one has the Schwarzschild/non- 
Schwarzschild black-hole solutions of Section |9.3| plus the ‘wormhole’ solutions of Section 


9.2 When the latter are Z2 symmetric, they lead cleanly through to a second sheet of 


spacetime, as in Eq. (9.2). The asymptotically flat (2,2) solutions, both Schwarzschild and 


non-Schwarzschild black holes, and also the Z2 wormhole solution all appear as separatrices 
between more generic singular solutions found numerically. The nature and interpretation 
of the latter remain to be better understood. 


10.1 Stability issues for black holes 


Going beyond the static spherically-symmetric ansatz (1.3) is outside the scope of the 


present paper. But one can contemplate what could happen dynamically once time depen¬ 
dence is allowed. A full stability analysis of the various phases of the static solution space 
would be desired, but in the meantime one can extract some partial stability information 
from various quasinormal mode studies of the stability of the Schwarzschild solution itself, 


considered as a solution of the higher-derivative (1.1) theory. This has been studied, e.g., 
in Ref. [22]. It was found there that the Schwarzschild solution is stable in the 7 R + (3R 2 
theory with a = 0. This is not surprising, because that theory is equivalent o [23] to 
ordinary general relativity coupled to a positive-energy massive scalar field. 

In Ref. [22] it was also suggested that the Schwarzschild solution could become unstable, 
for nontachyonic values of (m 2) 2 = 57, for sufficiently small values of 


g,w — 


Mm ,2 
Mp\ 


( 10 . 1 ) 


where Mp\ is the Planck mass. Ref. |22j then went on to claim, nonetheless, that detailed 


analysis of the quasinormal modes of the theory (1.1) showed no such instability. This 
conclusion has, however, been challenged more recently in Ref. [24] . where it is claimed 
that Ref. [2)2] erred in considering only a static S-wave potential instability. Instead, the 
analysis of Ref. [24] does fold Schwarzschild S-wave instabilities for < 1 by treating the 
Ricci tensor R pv as an effective massive field. This instability is compared to Schwarzschild 


instabilities found in massive theories of gravity [25] . 

Instability of the Schwarzschild solution for small black holes (i.e. small g, w ) raises 
the question whether a stable sector of the static solution space exists, and whether one 
or another of the non-Schwarzschild solutions we have discussed could then represent a 
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stable final phase. Clearly, the relation between and the branch point in the black- 
hole solution space could be an important issue in this regard. The analysis of time- 
dependent gravitational collapse can, however, be quite involved, as indeed it is already for 
the apparently simpler case of Brans-Dicke theory [26]. 
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Appendices 


A Reduced nonlinear field equations 


The reduced field equations of maximal third order (2.14) as derived in Section 2.1 for 
spherically-symmetric solutions written in Schwarzschild coordinates are given in detail as 
follows, where A^ and B^ are third derivatives. For equation (2.14a) we have 


24r 4 A 3 B 4 H rr = 8r 3 A 2 B 2 B^ ( r(a - 3/3)5' - 2(a + 60)B) 

- 4r 2 AB 2 A" ( r 2 (a - 3/3)5' 2 - 4r(a + 6/3)55' + 4 (a - 12/3 )B 2 ) 

- 4r 4 (a-3f3)A 2 B 2 B" 2 

- 4r 2 ABB" ^2 rBA! (r(a - 3/3)5' - 2(a + 6/3)5) 

+ A (3 r 2 (a - 3/3)5' 2 - 12r(a + 3/3)55' + 8(a + 6/3)5 2 ) ^ 

+ 7 r 2 B 2 A' 2 ( r 2 (a - 3/3)5' 2 - 4r(a + 6/3)55' + 4(a - 12/3)5 2 ) 

+ 2r 2 ABA'B' (3r 2 (a - 3/3)5' 2 - 4r(2a + 3/3)55' + 4(a + 24/3)5 2 ) 

+ 24A 3 5 3 ( 7 r 3 5' + 5 ( 7 r 2 - 12/3)) 

+ A 2 (V 4 (a - 3/3)5' 4 - 4r 3 (5a + 12/3)55' 3 

- 4r 2 (a - 48/3)5 2 5' 2 + 32r(a + 6/3)5 3 5' - 16(a - 21/3)5 4 ^ 

+ 8vl 4 5 4 (2a - 6/3 - 3 7 r 2 ) , (A.l) 
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while for equation (2.14b), using the definitions of X(r ) and Y(r) as given in (2.13), we 
have 


2 r 4 A 5 B 2 ( arB ' - 30r5' - 2 aB - 12/3B ) 2 [H tt - X(r)H rr - Y(r)d r H rr ) = 

72a0r 3 A 2 A (3) 5 4 (r(a - 3/3)5' - 2(a + 6/3)5) 

+ 3Qa(3r 2 AB 3 A" ^13 rBA! (2/a + 6/3)5 - r(a - 3/3)5') 

- 2A(-r 2 (a - 30)5 ' 2 + r(a + 6/3)55' + 2(a + 60)5 2 )^ 

+ 120r 4 (a - 3 f3)A 3 B 2 B" 2 ((a + 6/3)5 - r(a - 30)5') 

+ Ar 3 A 2 BB" ^3(3BA' (r 2 (a - 30) 2 5 ' 2 + r (a 2 - 15a0 + 36/3 2 ) BB' - 6 a(a + 6/3)5 2 ) 

- 30AB' (—r 2 (a - 30) 2 5 ' 2 - 6 ar(a - 3/3)55' + 2 (7a 2 + 48a/3 + 36/3 2 ) 5 2 ) 

+ 7 (—r)(a - 3(3)A 2 B 2 (2(a + 6/3)5 - r(a - 3/3)5') ^ 

+ 5O4a0r 3 5 4 A ' 3 (r(a - 3/3)5' - 2(a + 6/3)5) 

- 3/3r 2 AB 2 A' 2 (V 3 (a - 3/3) 2 5 ' 3 + 3r 2 (l7a 2 - 57a/3 + 18/3 2 ) 55 ' 2 

- 60a?’(a + 6/3)5 2 5' - 4 (23a 2 + 150a/3 + 72/3 2 ) 5 3 ^) 

- 60?n4 2 5A'(V 4 (a - 3/3) 2 5 ' 4 + r 3 (lla 2 - 39a/3 + 18/3 2 ) 55 ' 3 - 4r 2 ( 8 a 2 + 51a/3 + 18/3 2 ) 5 2 5 ' 2 
+ 4r (lla 2 - 12a/3 + 18/3 2 ) 5 3 5' - 16 (4a 2 + 21a/3 - 18/3 2 ) 5 4 ^) 

+ A 3 ^ - 4r(a - 3/3)5 4 5' (l2/3(5a + 3/3) + r(a - 3/3)A' ( 7 r 2 - 12/3)) 

- 2r 2 5 3 5 ' 2 (6/3 (a 2 + 66a/3 + 36/3 2 ) + 7 r 3 (a - 30) 2 A') 

- 8 (a + 6/3)5 5 (—6/3(5a + 3/3) — rA! ( 2 a ( 7 r 2 — 6/3) + 3/3 (120 + 7 ^ 2 ))) 

- 30r 5 (a - 3/3) 2 5 ' 5 + 30r 4 (-19a 2 + 51a0 + 180 2 ) 55 ' 4 + 120r 3 (l3a 2 + 84a0 + 36/3 2 ) 5 2 5 ' 3 

- 8 A 5 B 4 ^ r(a — 30)5' (a ( 7 r 2 — 60) + 60 (30 + 7 r 2 )) 

+ (a + 60)5 (a (60 - 27 r 2 ) - 3(3 (60 + 7 r 2 )) ^ 

- 2A 4 B 2 ^ 7 r 5 (a - 30) 2 5 ' 3 - 6 r 2 (a - 30)55 ' 2 (a ( 7 r 2 - 40) + 30 (40 + 7 r 2 )) 

+ 4r(a - 30)5 2 5' (a ( 7 r 2 - 240) + 60 ( 7 r 2 - 60)) + 4 (2a 2 + 15a0 + 180 2 ) 5 3 (120 + 7 r 2 ) ^ . 

(A.2) 
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From these two coupled third-order differential equations, one anticipates that the so¬ 
lution will depend in general on a total of six integration constants. For a pair of linear 
differential equations, this can be demonstrated straightforwardly by reducing the system 
to a single sixth-order differential equation for just one function, e.g. A(r), by repeatedly 
differentiating and substituting between equations so as to eliminate B{r ) and its deriva¬ 


tives. In the present highly nonlinear equation system (A. 1, A.2), this is not feasible to do 


explicitly because this would involve the inversion of polynomials of order > 4. However, 


the idea can be outlined as a sequence of operations on Equations (A. 1, A.2) as follows: 


(Al) :0 = 


(A.2) : 0 = gi(r, A, B, A', B’, A', B", A'") 


6 > T .(ph 2 |) : 0 = d rgi (r, A, B, A, B’, A", B", A") 

= g 2 (r, A, B, A, B’, A', B", A", A", A 4 )) 
.-. B"' = 9 2 1 (r. A, B, A, B\ A', B", A", A< 4 >) 


sub into (Al[ : 0 = / 2 (r, A, B, A, B', A', B\ A", A A) 
■ B" = / 2 _ 1 (r, A, B, A, B', A", A", A 4) ) 


sub into (A.2) :0 = gs(r, A, B, A, B', A”, A'", A 4 )) 
.-. A= 5 3 1 (r,A,B,A,A , ,A' , ,A( 4 )) 


sub / 2 1 and 53 1 into (A.2) : 0 = 54 (r, A, B, A, A', A", A 4 ^) 

A^^r.V'AU"',#) 


sub into (A.2) : 0 = g 5 (r, A, A, A', A", A (4 \ A®, A< 6 >). 


(A.3) 


B Coupling of a shell source to the higher-derivative theory 

B.l Coupling an (0,0) vacuum inside the shell to a (2,2) solution outside 

One can carry out a successful coupling of a thin-shell delta-function stress-tensor source 
to the full nonlinear higher-derivative theory in a fashion similar to the couplings in the 
linearised theory as discussed in Section [4j As we saw in Section [5j it’s only with an 
exterior (2, 2) family solution combined with an interior (0,0) family vacuum solution that 
the count of available interior plus exterior solution parameters is sufficient to satisfy the 
nine continuity, step, asymptotic flatness and asymptotic Minkowskian requirements. In 
this appendix, we discuss in more detail how these coupling requirements can be met. 

Coupling a thin shell source to an external (2,2) family solution faces a number of 
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challenges. Principal among these is the lack of a closed-form expression for the (2, 2) 
family of solutions, so one is limited to carrying out the coupling using series solutions as 
given in Section [3j There is also an awkwardness arising from the choice of the Schwarzschild 
form ( |1.3[ ) for the metric. Such difficulties were already noted in the classic treatment of 
delta-function couplings in general relativity given in Ref. m- Such difficulties might be 
alleviated by working in other than Schwarzschild coordinates, but we have not explored 
this possibility in detail. The difficulty with the in «->• out matching in the series-expanded 
theory is a leading-term matching problem in A: the inner solution has j 4(0) = 1 and so we 
would naively expect A(£_) « 1 + 0(£ 2 ) but the outer solution has A{£ + ) « c ^£ 2 + 0(£ 3 ), 
while A is continuous so A(£_) and A{£+) need to be equal for arbitrary £. The £ —> 0 
limit is particularly interesting because this will give us the field of a point source in the 
higher-derivative theory. In order to make this work for arbitrary small £ we need to elevate 
the free parameters of the solution to functions of £. 

Let us first illustrate the method with general relativity, where such dependence on the 
shell size is well-known [12], 


B.2 Shell coupling in general relativity 

To set the stage for the more involved coupling problem in the higher-derivative theory, we 
first review the analogous coupling problem for distributional sources in general relativity. 
Some classic references for this are [T2] and [13]. 


B.2.1 Shell-source coupling using the closed-form Schwarzschild solution 


Consider general relativity with a thin-shell source as in (5.1 5.3), with T tt = ^5(r-£). 
It is convenient to define the length scale Lm '■= 2 GM = M(8iT'y)^ 1 . Coupling the source 
to the equations of motion, one can show that B is continuous while A has a step: 

Lm ^m (£-) 2 


B out (£ + ) = B in (£_) , 


A out (£ + ) - A in (£_) = 


£B(£)-L M A in (£_) ' 
In terms of these parameters the held of a spherical shell is 


(B.l) 


— 1 ; 
Bin = b , 


Aant, — 


1 - w 1 

br 


Be, lit — 


1 _ Lm 
1 be 


Lm \ 

br ) 


(B.2a) 

(B.2b) 

(B.2c) 


(B.2d) 
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where the exterior solution has the form of the Schwarzschild solution 

^4out = ~ yy , (B.3a) 

r 

B out =k 2 , (B.3b) 

provided the interior free parameter b scales and the Schwarzschild radius r s is related to 
the source length L\f via the f-dependent expressions 

b=k 2 [ l-j), (B.4a) 

L M =k 2 (l- 1 ^)r s . (B.4b) 

Before proceeding with the core of our discussion, the signs here need a comment. For the 
case 0 < £ < r s a horizon exists and M and b are negative. It is familiar fact that we 

have —b ++ signature outside the Schwarzschild horizon (r > r s ), and signature H-b+ 

inside the Schwarzschild horizon (£ < r < r s ). At the shell source (r = £), the function B 
is continuous and non-vanishing so it is therefore positive inside and out. The function A, 
however, has a step at r = £. Outside the source (for £ < r < r s ), A is of course negative, but 
inside the shell source (for r < £) the equations of motion in the Schwarzschild-coordinate 
metric ansatz (where the angular part of the ds 2 metric is just r 2 dkl 2 ) require A = 1. In 
consequence, the signature inside the shell is + + ++. The fact that the source mass M 
has opposite sign to that of the Schwarzschild radius r s is not unexpected - the source is 
static in a region where t is spacelike, i.e. the source is of tachyonic sign. In the higher- 
derivative case this peculiarity will not arise, because we shall find that there is no horizon 
in the source-coupled solution, and the metric components A and B are continuous across 
the shell source. We thus anticipate a —b ++ signature for all r in the higher-derivative 
theory. 

A key point in the above analysis is the fact that in general relativity the interior free 
parameter b has to blow up as l~ l as the shell is shrunk down to a point. We shall find that 
the higher-derivative case also requires poles at £ = 0 in the free parameters of the interior 
solution. 

B.2.2 The Schwarzschild solution from a series-solution point of view 

To set the scene for analysis of coupling in the higher-derivative theory, where series solutions 
will be all that we have available, let us now repeat the above coupling discussion for the 
Einstein-theory Schwarzschild solution using only a series solution. The interior series 


49 



vacuum solution is of (0, 0) structure: 

^4(0,0) = 1 + ... (B.5a) 

B (o,o) = b + ... (B.5b) 

and the solution exterior to the source is of (1,-1) structure: 

= xr — x 2 r 2 + x 3 ?’ 3 — x 4 r 4 + 0(r 5 ) (B.6a) 

= ^+xy + .... (B.6b) 


To implement the matching conditions, we need to elevate the free parameters (x, y, b ) to 
functions of £: (x(£),y(£),b(£)). 

In order for the solution exterior to the source to be of the (1,-1) family, we let the 
free parameters be expressed as Taylor series: 

x{£) = x(0) + £x'(0) + ^£ 2 x"(0) + ... (B.7a) 

y(£) = 3/(0) + £y\ 0) + \py"( 0) + ... . (B.7b) 

For the interior free parameter b, however, we need to use a Laurent series - i.e. we allow 
£~ n poles: 

b = £ a (bo + hi + M 2 + b 3 e 3 + ...) . (B.8) 

Then, in order for y(£) to remain finite in the l —> 0 limit, we find that we need to scale 
Lm (Le. we scale the mass M) as well: 

L m = £ d (T 0 + £L\ + £ 2 L 2 + £ 3 L 3 + ...). (B.9) 

We find that the suitable poles have a = — 1, d = — 1. 

The solution is: 

y(£) = - L m (£) {£x(£) - £ 2 x (£) 2 + O (£ 3 )) 

= - L 0 x( 0) + £ (T 0 (x(0) 2 - x'(0)) - Tix(0)) + O (£ 2 ) (B.lOa) 

b{£)= - L m (£)x(£) + O (£ 3 ) 

= ~ L ° 1 ^ > + (-Tox'(O) - Tix(0)) 

+ h (-T 0 x"(0) - 2Tix'(0) - 2T 2 x(0)) + O (£ 2 ) , (B.lOb) 

which matches the analytic solution (where we hold r s fixed, something we don’t know how 
to do in a series solution) for small £ upon renaming x(0) = — y- and Lm( 0) = —k 2 r 2 , 
and which also can be solved by matching for all £, producing x(£) = —T and Lm {£) = 
k 2 (1 - r f) r s . 
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B.3 Shell-source coupling in the higher-derivative theory 


Now we consider the coupling of a thin-shell stress-tensor source of the form (5.1 - 5.3) to the 


higher-derivative theory (1.1). Inside the shell source, we require a (0,0) vacuum solution, 
and outside the shell we consider a ( 2 , 2 ) solution with the following notation: 
r 3 v 3 w 2 r A (u > 2 (2v 2 {v 2 W 2 - 4u 4 ) + v 2 )) 


A 2 i u i w 2 • y^zy^^zy^z^z i ^3 J 5 

A = r W 2 H--- —9 -— + r w 5 + 0(r ), 


V2 

B = r 2 v 2 + r 3 v 3 + r 4 n 4 + r 5 v 5 -f 0(r 6 ). 

The form (for at least n < 14) of the (0,0) solution is 

A =1 + 02 T 2 + 'y ' 


6r 2 


r" I Z 


n,p,q,m 


P 




a , 


B 

Vo 


: 1 + hr 2 + Y n , P ,m, q r n ( p 


v-P 


<«■ 


(B.lla) 

(B.llb) 

(B.12a) 

(B.12b) 


n,p,q,m 

where the X niP _ m _ q and Y n , Ptrrhq are rational numbers and the n,p, q, m sums are taken over 
n = 4, 6 , 8 ,...; 1 < p < 7 ^', 0 < q < ^ — 1 and 0 < m < p. 


Similarly to the method used for general relativity, (B.7a), (B.7b), (B. 8 ) and (B.9), we 
elevate the free parameters to functions of £ in the following scheme 


b 0 =£ 2 H 0 (£), 

W 2 =W 2 (£) , 

a2=r 2 G 2 (£), 

V2 =V 2 (£) , 

b 2 =£~ 2 F 2 (£), 

V3 =V 3 (£) , 


II 


u>5 =W 5 (£) , 

M=£ d p {£), 

V5 =V 5 (£) , 


(B.13) 


where we have factored out the poles so the remaining functions are Taylor series in £. 

To carry out matching across the source shell, we need to have A- m (£.) = A out (£ + ), 
where 

\k / 0 \ q 

, (B.14) 


A iD (£_) = l + G 2 (£)+ < 

k,n,q,m 


7 


P 


n — k 


,m,q 


G2{£) m F 2 {£)^- m ^ 


a 


where the sum is taken over k = 0, 2,4, 6 , ...; k + 2 < n = 4, 6 , 8,...; 0 < q < f — 1 
and 0 < m < Gathering powers of £, we have An(^-) ~ £°+ ~ £ lj r ~ £ 2 + ... and 

A out (£ + ) ~ £ 2 + ..., so the leading-term matching problem in A is now displayed in the 
vanishing of the interior £° term, where one has for Aq (the £° term in An) 

A) = 1 + G 2 (0) + X nthm>q G 2 (0) m F2(0)5- m (£)' , (B.15) 

n,q,m ' ' 
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where the sum is taken for n = 4,6, 8, ■■■', 0 < q < S — 1 and 0 < rri < S. The higher terms 


in (B.15) need to decrease in amplitude with n so that the sum converges, and its limiting 
value as n -+ oo must vanish so as to match the structure of A out . We have not carried 
out an exhaustive analysis of the convergence properties of the resulting series, but we may 
consider the structure in the simplifying limit /3 « a. In this limit, the q > 1 terms are 
suppressed and we need only consider X n ^ m 0 . In this limit, the numbers X n n n 0 are 
equal to 1 for all n while the other X n n m< n_ l 0 grow at most linearly with n (for at least 
n < 14). If we rename t = ^n, X n ^ m 0 = X t>rn and G 2 {0) = £F 2 (0) and take the sum only 
out to T terms, we have 


An = 


Xt,mC m F 2 ( Of- 


(B.16) 


*=1,2,3,4,...T 

0 <m<t 


Given the at most linear growth in t of Xj >m , one gets an estimate of the sum Aq by letting 
Xt t m = a + bt, for which an estimate sum Aq can be carried out: 


A 0 = 


E 


1 - C t+1 

0 a + bt ) - t ^-F 2 (0)* ! 


(B.17) 


t=l,2,3,4,„.T 

for which the ratio of successive terms at large t is 

a + b + bt F 2 ( 0) t+1 1-C t+2 
a + bt F 2 (0y 1 - C t+1 


f 2 (o) ■ 


(B.18) 


Accordingly, convergence is obtained in the (3 « a limit. Convergence is expected when 
the f3 « a limit is relaxed as well. 

Once the 1° terms are matched inside and out, the matching of l N - 1 terms and the 
matching of A'(£), B (£), B'(£), B" (£), A" (£) should follow suit with less difficulty. The final 
result will describe the interior and exterior metrics in terms of ct, /3,7, Lm{£) and three 
other free parameters Pi{£),p- 2 (£),P 3 {£)■ 


C On higher-derivative no-hair theorems 

In this appendix, we present a recalculation of the argument of Ref. |14| including a cos¬ 
mological constant. Unfortunately, this corrects the calculation of [Til in a way that voids 
its conclusion about the no-hair consequences of the non-trace part of the higher-derivative 
gravity field equation. 

The discussion of [llj can be generalised to the case with a cosmological constant. Take 
as Lagrangian 

I = J d A x^g ( 7 (i? - 2A) - aC lwprJ C lwp(T + PR 2 ) . (C.l) 
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The equation of motion (2.1) gains a term 


making the trace 


where we have defined 


H>; = 6/3 (HR — ml(R — 4A)) 

= 6/3 (ns - mis ) , 

S = R- 4A. 


(C.2) 

(C.3) 

(C.4) 

(C.5) 


For the trace no-hair theorem, the discussion proceeds from this point on just as in the case 
without a cosmological constant as reviewed in Section [6j obtaining finally 

0 = J Vh d 3 x [D a (A S D a S) - A ( D a S ) ( D a S) - A m 0 2 S 2 ] . (C.6) 

The outer boundary contribution vanishes if D a S = D a R falls off appropriately fast at 
infinity. Consequently, we deduce that 


5 = 0 4^ R = 4A (C.7) 

if the inner boundary term vanishes, which is ensured if the inner boundary is at a horizon. 

The argument for the traceless part of the higher-derivative no-hair theorem runs into 
trouble, however. 

Define the shifted quantities 


S jW = R, lv - g llv A, S = <rS lw = R- 4A. 


(C.8) 


Letting denote the Ricci scalar of the spatial part of the metric h a b , we also define the 
shifted quantity 

(3)5 = ^R - 2A 


and finally we define 7712 (A) such that 777 - 2 ( 0 ) = m 2 

7772 (A ) 2 : = 


2 _ 7 + |A(3/3 — a) 


2 a 


thus obtaining the equations of motion 

H, 


0 = 


L His 


—2a 


(C.9) 


(C.10) 


5=0 


= □ V - 777 2 (A) 2 5 M j, + 2 S?Sy P - 2V P V M 5£ - - 9l j,yS pa Sp a , (C.ll) 
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which smoothly analytically continue from A = 0 to A finite. Note that the contracted 
Bianchi identity in this case is 


= - V U R , 


V /i 5 #H/ =0 for 5 = 0. 


1 ( 3 )5 

V^5 M4 = 0 = LfiSij + -SijD3 A + — A A. 


(C.12) 

(C.13) 

(C.14) 


At this point, we can contract (C.ll) with A 5^", dimensionally reduce on the t coordinate, 


and integrate over three-dimensional space while using (C.7) to obtain 


0 = / Vh d 3 x 

A 5^ 


Js 

—2a 

5=0. 


(C.15a) 


= j v7i d 3 x A Q ( 3 A A f 3 )^ + A S kl D i S kl - 2A S k iD k S li - A ^SDjS j *\ (C.15b) 

(C.15c) 


A 


- - A A 5 A (3) 5 + 2A A (3 )SD j S ji - A AA fe [A 5^ - 2A'5fci] 


- A 


( 3 ) 5 ' 


(m 2 (A) 2 + (3)5) - A S ij Sij (m 2 (A) 2 - 25) 


(C.15d) 


Unfortunately, owing to the fact that the squared terms in the middle line (C.15c) of the 


integrand are of non-uniform sign, it is not possible in the case of higher-derivative gravity 
to say that the assumption of a horizon and suitable boundary conditions at infinity imply 


that Sfu, vanishes, even for large m 2 (A). This continues smoothly to A = 0, where (C.15) 
with a horizon and asymptotic flatness unfortunately does not imply that vanishes, 
even for large m 2 . 
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